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ABSTRACT 

In this thesis, the theory of associated polynomials is 
applied to linear circuits and systems. Associated polyno- 
mials are defined and used in the presentation of a computer 
algorithm suitable for accomplishing a partial fraction ex- 
pansion including repeated roots. Variations of the coeffi- 
cients of a polynomial are related to the variation of its 
roots (and vice-versa) by associated polynomials. These re- 
sults are used in a root solving process and to exoress sen- 
sitivity coefficients in; a new analytical form. Lastly, a 
theory of compensating parameter adjustments is developed 
with applications to self-adaptive systems. It is shown that 
by compensating adjustments of any two parameters in the sys- 
tem which are linearly related to the coefficients of the 
characteristic equation, adaptive compensation is possible 
for a real root of the closed loop transfer function only 
and that for complex roots, the real oY imaginary parts may 
be kept invarient but not both. 



fJRY 

VAL POST 
ON'CERF.Y , 

Chapter 

I 



II 



III 



IV 

V 



GKADUA 

CALIF. 



L'fi .iChuJL 
93940 



TABLE OF CONTENTS 
Section Title 

Associated Polynomial Theory 

1.1 Introduction 

1.2 Definition of Associated Polynomials 

1.3 Characteristics of Associated 

Polynomials 

Partial Fraction Expansion Using 
Associated Polynomials 

2.1 Introduction 

2.2 Case (1): Partial Fraction Expansion 

2.3 Case (2); Partial Fraction Expansion 

2.4 Case (3); Partial Fraction Expansion 

Polynomial Sensitivity Theory 

3.1 Introduction 

3.2 Coefficient Sensitivity 

3.3 The Basic Orthonormality Relation 

3.4 Root Sensitivity 

3.5 Root Solving Using Associated 

Polynomials 

Sensitivity Coefficients 

4 . 1 Introduction 

4.2 Theory of Sensitivity Coefficients 

A Two Parameter Adaptive System Using 
Associated Polynomials 

5.1 Introduction 

5.2 Basic Theory 

5.3 Case (1) ; The Dominant Root is Real 

and is to be Held Constant 



Page 

9 

9 

11 

15 

23 

23 

24 

27 

28 
38 
38 
38 
43 
49 

51 

56 

56 

57 

62 

62 

66 

69 



3 



Chapter Section Title Page 

5.4 Case (2) : The Dominant Roots are a 

Complex Conjugate Pair and are to 

be Held Constant 75 

5.5 Case (3) : The Dominant Roots are a 

Complex Conjugate Pair and the 
Damping Coefficient is to be Held 
Constant 83 

5.6 Case (4): The Dominant Roots are a 

Complex Conjugate Pair and the 

Damped Resonant Frequency is to be 

Held Constant 87 

VI Conclusions 91 

Bibliography 93 

APPENDIX A Horner's Method 95 

A.l Introduction 95 

A. 2 The Nested Evaluation of Polynomials 95 

A. 3 Synthetic Division 96 

A. 4 Approximation of the Roots of a 

Polynomial Equation 98 



4 



LIST OF ILLUSTRATIONS 

Figure Page 

4.1 Dynamic System Basic Block Diagram 57 

5.1 Feedback Control System 62 

5.2 Control System with variable parameter 6 64 

5.3 Self-Adaptive Control System 65 

5.4 Block Diagram Example 5-1 System 70 

5.5 System Poles With a Variable Parameter 71 

5.6 Block Diagram Example 5-2 System 74 

5.7 Change in Root Location With a Change 

in a System Parameter 76 

— +• H 

5.8 n Order Feedback Control System 79 

5.9 Modified Closed-loop Feedback Control 

System 80 

5.10 (n-1) Order Feedback Control System 82 



5 



LIST OF TABLES 



Table 

5.1 



Page 

System Roots With a Variable Parameter 70 



7 



ACKNOWLEDGMENT 



I would like to express my appreciation to Professor 
S. R. Parker, for the encouragement and assistance in the 
preparation of this thesis. 



8 



CHAPTER I 



Associated Polynomial Theory 
1.1; Introduction 

This thesis is concerned with the investigation and 
application of associated polynomials to the analysis and 
synthesis of linear systems. As is shown, every polyno- 
mial of degree n can- be used to generate a sequence of n-^l 
polynomials of decreasing degree called the associated 
polynomials of the defining polynomial. The associated 
polynomials depend in a particularly simple way upon the 
coefficients of the original polynomial. Due to the rela- 
tionship between the coefficients of the original polyno- 
mial, the associated polynomials, and the properties of the 
original polynomial which can be expressed simply in terms 
of the associated polynomials, several useful aoplications 
and results are possible in network analysis and control 
theory. 

The concept and use of associated polynomials is not new 

in that these polynomials have appeared in many guises in the 

literature , usually related to some particular application. 

Associated polynomials first appeared in the literature in 

the late eighteenth century when Legrange used them to obtain 

( 2 ) 

a solution of the homogeneous linear difference equation 

(1 3) 

Associated polynomials have been called Horner polynomials ' 
as they were used directly in Horner's method or rule^^^. 
the purely mathematical aspects of associated polynomials 
not of direct significance to this thesis. The use of asso- 
ciated polynomials in Horner's method is presented separately 
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in appendix A. What often occurs in engineering applications 
is not an associated polynomial itself, but rather a number 
corresponding to an associated polynomial evaluated at a 
particular root of the original polynomial. 

Recently^^^, the properties of associated polynomials 
have been investigated in themselves and new concepts and re- 
sults introduced. Reference 1 is written in a pure mathemat- 
ical language with only mathematical formulas as results. 

The object of this thesis is to demonstrate that associated 
polynomials have very useful properties and applications in 
the fields of circuit and systems theory. 

This thesis, divided into five chapters, presents the 
basic associated polynomial theory in chapter one. In chapter 
two, this basic theory is used to present an original algo- 
rithm suitable for computer programming for accomplishing a 
partial fraction expansion including repeated roots of high 
order of multiplicity. Two theorems are developed in chapter 
three that show that the variation of the coefficients of a 
polynomial with the variation of its roots or vice-versa is 
determined by associated polynomials. The theorems are made 
possible by deriving the orthnormality relationship between 
the set of associated polynomials and the roots of the defin- 
ing polynomial. The above theorems have important engineering 
applications in that they are used in a root solving process 
in chapter three, and in chapter four to express sensitivity 
coefficients in a new analytical form. Chapter five shows 
how associated polynomial theory may be used in the develop- 
ment of a theory of compensating parameter adjustments with 
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application to self-adaptive systems in which a change in a 
variable parameter in a plant may be compensated for bv a 
change in another variable parameter such that the dominant 
root performance of the overall system is controlled. An 
important theorem of this chapter is the proof that if a feed- 
back control system has a real root then the location of this 
root may be kept constant by compensating adjustments of any 
two parameters in the system which are linearly related to 
the coefficients of the characteristic equation. It is shown 
that when the compensating parameters are associated with 
simple zeros and poles of the open loop transfer function, 
such adaptive compensation is possible for a real root of the 
closed loop transfer function only and that for complex con- 
jugate roots the real or imaginary parts may be kept invar- 
ient but not both. This theorem is of fundamental importance 
in adaptive control theory. 

1.2; Definition of Associated Polynomials 

Let a monic polynomial, F(s) , of degree n be expressed 

as ; 



The condition that F(s) be monic, i.e., a^ = 1, is convenient 
rather than essential. Let the distinct roots of F(s) be 
labeled, p ^ , j = 1,2,. ..,n. 

If u is any variable different from s, then the oolyno- 



F(s) 



a s 
n 



n 




1 ( 1 - 1 ) 



mials Aj (u) ; j = 0,1 



/ • • • / 



n are defined by the operation: 



sF (s) -uF (u) 
s-u 



n 



I A. (u) s^ 
= n ^ 



= Aq (u) + Aj^(u)s + ... + A^(u)s^ 



( 1 - 2 ) 
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The set of polynomials, (u) ; j = 0,1,..., n are called the 
associated polynomials of F(s) . 

Example 1-1 

2 

Given the polynomial: F(s)=s +3s+2 

The associated polynomials of F(s) are found by the following 
operations : 

F(u) = u^ + 3u + 2 

Performing the multiplication and division of equation 1-2, 

sF(s)-uF(u) 2 , . , , 2,^ 

^ = s + (3+u)s + (u +3u+2) 

s-u 

The associated polynomials of F(s) are: 

Aq(u) = u^ + 3u + 2 
Aj^(u) = 3 + u 
A2(u) = 1 

The basic definition of associated polynomials given by 
equation 1-2 may be simplified by noting the relationship of 
associated polynomials and the coefficients of the basic 
polynomial. This relationship is given by Theorem 1-1. 

THEOREM 1-1 

th 

The associated polynomials of an n order polynomial, 
F(s) , may be foxmd from the coefficients of F(s) by the recur- 
sion formula: 

n . _ . 

A.(u) = ^ a.u^ ^ : j = 0,1,. ..,n 

^ i=j ^ 
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Proof: 



Multiplying equation 1-2 by (s-u) ; 



sF(s)-uF(u) = [A^ (u) +A, (u) s+. . -+A (u)s’^](s-u) (1-3) 

u 1 n 

Substituting equation 1-1 for F(s) and multiplying through 
by (s-u), equation 1-3 becomes: 

as + a ,s + ... + a,.s - uF(u) = A (u)s + 
n n-1 u n 

[A^_j^ (u) -uAj^(u) ] s^+ + [Aq (u) -uAj^ (u) ]s - uAq(u) 

(1-4) 



Comparing the coefficients of the equal powers of s, 



^n 


= A (u) 
n 




(l-5a) 




= A , (u) 


- uA (u) 


(l-5b) 


n-1 


n-1 


n 


a o 
n-2 


= '^n-2 <“> 


- uA„ , (u) 
n— 1 


(l-5c) 



^0 ~ “ uAj^(u) (l-5d) 

Since a^^ = 1, equation 1-5 may be explicitly expressed as: 

Aj (u) = uAj^j^(u) + S’j 5 j = 0,l,...,n-l (1-6) 

By definition of the set of associated polynomials, A. (u), j 
ranges from 0 to n. Hence let: 



Using equation 1-7, 



Vi<“> = 



equation 1-6 



0 

becomes : 



(1-7) 



Aj (u) = uAj^jj^(u) + : j = 0,1,. ..,n 



( 1 - 8 ) 
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In order to form an explicit expression for the associated 
polynomials, note: 



A (u) 
n 




(l-9a) 


A , (u) 
n-1 


= uA(u)-ha T=a ,+au 
n n-1 n-1 n 


(l-9b) 


A , ( u) 
n- 2 


.V 2 

= uA^ T (u) + a^ = a^ + a^ tU + a^u 

n-1 n-2 n-2 n-1 n 


(l-9c) 


Aq (u) 


= uA, (u) -f a^ = + a, u + ... + a 

1 U U 1 n 


(l-9n) 


for j = 0 : 




(1-lOa) 


for j = 1: 


A, (u) = 1 a.u^"^ 

^ i=l ^ 


(1-lOb) 


hence : 


n . _ . 

A.(u) = J a.u^ ^ j = 0,1,. ..,n 

^ i=j ^ 


(1-11) 



Corollary 1-1 

The associated polynomial, (u) , is a polynomial of 
degree n - j . 

Corollary 1-2 

The particular associated polynomial, Aq(u), is always: 

Aq(u) = F(u) (1-12) 



Example 1-2 

Using the results of Theorem 1-1, the associated polyno- 
mials of the polynomial, F(s) , of example 1-1 may be deter- 
mined as follows : 
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F(s) = s + 3s + 2 = 



^ 2 ® 



+ + Hq 



A. (u) = I : j - 0,1,2 

i=j 

The associated polynomials of F(s) are: 

A2(u) = a2 = 1 

A^(u) = a^^ + a 2 U = 3 + u 

2 2 

Ag(u) = ag + a^^u + a 2 U = 2 + 3u + u 
which agrees with example 1-1. 

1.3: Characteristics of Associated Polynomials 

The above section has shown the elementary relationship 
of associated polynomials to the coefficients of the original 
polynomial. There are other interesting properties of asso- 
ciated polynomials that can now be considered. These prop- 
erties include the relationship between associated oolyno- 
mials and the original polynomial which leads to the concept 
of synthetic division and the fact that associated polyno- 
mials are elementary symmetric functions. These properties ' 
are stated in the following theorems. 

Theorem 1-2 

The relationship between a polynomial, F(s) , and its 
associated polynomials, given any particular root, p^, of 
F ( s ) i s : 



F(s) 

s-p. 



I A . (p. ) S^ ^ 
j=l => ^ 
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Proof : 



For convenience, equation 1-2 is repeated as: 



(sF(s)) [s,u] = Aq(u) + Aj^(u)s + ... + A^(u)s^ (1-2) 



n-1. 



Equation 1-2 may be rewritten as: 

(sF (s) ) [s ,u] = Aq (u) + s [Aj^ (u) + A 2 (u) s + ... + (u) 

(1-13) 

However the expression in the brackets of equation 1-13 may 
be equated to : 

^ = A, (u) + A^(u)s + ... + A (u)s^ ^ = J -A. (u)s^~^ 
s-u 1 2 n • 1 • 1 

3=1 

(1-14) 

Since F(u) = Aq(u) , equation 1-14 may be rewritten as: 



F(s) _ ■^0 



s-u 



u) n 

3=1 



.j-1 



(1-15) 



Note that equation 1-15 states that the quotient. and remainder 
after division by a monic linear polynomial may be calculated 
recursively. As shown in appendix A, the nature of this re- 
cursion leads to synthetic division. 

Let u = p^ , be a root of the polynomial, F(s) , then. 



since F(p^) = 0 



Aq (u) 



s-u 



= F(u) 
u=p^ s-u 



= = 0 



(1-16) 



u=p. s-p. 



and equation 1-15 may be rewritten as : 



?A.(p.)s 

-Pj 3 



3-1 



s-p 



= 1 3 



(1-17) 



1 j=l 

Theorem 1-2 may be numerically verified by the follow- 
ing example. 
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Example 1-3 



Given ; 



F(s) = s'^ + lOs^ + 35s^ + 50s + 24 = (s+1) (s+2) (s+3) (s+4) 
Let; 

Pi = -1; P2 = -2; = -3; = -4. 

The coefficients of F(s) are: 

a^ = 24; a^ = 50; a2 = 35; a^ = 10; a^ = 1 . 

The associated polynomials of F(s) are: 

A4(u) = a^ = 1 

A^Cu) = a^ + a^u = 10 + u 

A2(u) = 32 + a^u + a^u^ = 35 + lOu + 

2 3 2 3 

A^(u) = ^2^ ^ ^3^ ^ ^4^ = 50 + 35u + lOu + u 

Note; The particular trivial associated polynomial Aq(u) is 
no longer required and is not calculated in future examples. 
Consider the root, P 2 / of F(s) , then: 

Ai(p2) = 12; A2(P2) = 19; ^2(P2) = 8; A^(p2) = 1. 
Expanding equation 1-17 and evaluating: 

Iffi = A3_(p2) + A2(P2)S + A3(p2)s^ + A4(p2)s^ 



= 12 + 19s + 8s^ + s^ 



As a check: 

(s^+8s^+19s+12) (s+2) = s^ + lOs^ + 35s^ + 50s + 24 as given. 



17 



Corollary 1-3 



Theorem 1-2 developed a relationship between a polynomial, 
F(s), and its associated polynomials given a particular root 
of F(s) . This relationship may be extended to polynomials 
with multiple roots so that the following formulas are avail- 
able . 

CASE (1) ; TWO REPEATED ROOTS 




F(s) 



n-2 




n-3 





.(1-18) 



Expressed as an explicit relationship: 



(s-p^) h=0 j=h+2 




(1-19) 



CASE (2) : THREE REPEATED ROOTS 



F(s) 



A (p . ) S 
n 



n-3 




n-4 





+ [A^ (p^) +2p^A^ (p^) + . . .+ (n-2) (p^) ^^n^^i^^ 



( 1 - 20 ) 



Expressed as an explicit relationship 



(s-p^)^ h=l j=h+2 




j-(h+2)^ ^ j (1-21) 

D ■ 1 



18 



CASE ( 3) ; FOUR REPEATED ROOTS 






(s-p^) 






+ ... + [A^(pj^)+3p^Aj(p^) + + tl+. . .+ (n-3) 

( 1 - 22 ) 

Expressed as an explicit relationship. 



F(s) 



n-3 . , n 

= I I { 1 +, 



(s-p^) h=l j=h+3 

Example 1-4 
Given; 



.+ [i-(h+2) ] (p^) 



(1-23) 



F(s) = s^ + 13s^ + 67s^ + 171s^ + 216s + 108 = (s+2)^(s+3)^ 

Let; Pj^ = -2 ; P 2 = -3. 

The associated polynomials of F(s) are; 

Aj(u) = 1 

A^(u) = 13 + u 

A3(u) = 67 + 13u + 

A2(u) = 171 + 67u + 13u^ + 

Aj^(u) = 216 + 171u + 67u^ + 13u^ + u'* 

Evaluating the associated polynomials at ; 

Aj(Pl) = 1; A^(pj^) = 11; A^Cp^^) = 45; A2(Pj_) = 81. 
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Evaluating the associated polynomials at P 2 : 



= 1; (p^ = 10; A2(P2) = 37. 

Applying equation 1-18; 

^ ^ — 2 = A^(pj^)s^ + [A^ (pj^) +Pj^A^ (p^) ] s^ + [A^ (pj^)+Pj^A^ (p^) + 

(s-Pj_) 

2 2 3 

Aj(Pq^)]s + [A2(Pq^) '^4 ^5 ^ ^1 ^ ^ 

Evaluating ; 

■■ „ = Is^ + [ll-2]s^ + [45-22+4]s + [81-90+44-8] 

(s-pj^) 



or : 



F(s) 3 , Q 2 , . 0-7 f ■■ 5\3 

= s + 9s + 27s + 27 = (s+3) 



(s-pj^) 



Applying equation 1-21: 

— 3 = A^(p 2 )s^ + [A^ (p 2 )+ 2 p 2 Aj (P 2 ) ] s + [A^ (P 2 ) + 2 p 2 A^ (P 2 ) 
(s-p2> 



+3p2 A 3 (P 2 ) ] 



Evaluating : 



= Is^ + [10-6]s + [37-60+27] 



(S-P 2 ) 



or: 



+ 4s + 4 = (s+2) ^ 



(S-P 2 ) 



Theorem 1-3 

Any result which involves an associated polynomial, 

Aj (p^) , evaluated at a particular root, can always be inter- 
preted as a result involving elementary symmetric functions. 
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i . e . , 

■“ th 

where denotes the n-j elementary symmetric function of 

the roots of F(s) with absent. 

Proof : 

(5) 

Elementary symmetric functions are defined as fol- 
lows: If f(x) is a polynomial of degree n with complex coef- 

ficients, the functions obtained by taking the sum of all the 
roots of f(x) , the sum of all products of pairs of the roots 
of f(x) , the sum of all products of triplets of the roots of 
f(x) and so on to the product of all the roots of f(x) may 
be expressed rationally in terms of the coefficients of f(x). 
These functions of the roots of a polynomial, f (x) , are 
galled elementary symmetric functions. 

Hence, for a polynomial, F(s), of degree n with roots. 

Pi ' P2 > ' ' ' ' Pji • 

F(s) = (s-pj^) (S-P 2 ) . • • (s-p^) = s - a^s + 02 ^ -...+ (-1) 





(1-24) 


1 = Pi + P 2 + • • • + Pn 


(l-25a) 


2 ^ P 1 P 2 P 1 P 3 * * *. "^ Pn-l^n 


(l-25b) 


• • 

n ^ PlP2P3* * *Pn 


(l-25c) 



^1' ^2' ••*' elementary symmetric function of the 

roots , p^, P 2 , ... ,p^. 
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Rewriting equation 1-17 as: 



G (s) = (p . ) + A, (p. ) s + ... + A (p. ) ^ (1-26) 

s-p^ 1 ^1 2 n - 1 

G(s) is a polynomial of degree n-1. By comparison with equa- 
tion 1-24, the associated polynomials of F(s), evaluated at 
a particular root, p^ , can be seen to be elementary symmetric 
functions of p^, P 2 / . . . , p^ with the particular root, p^ = 0. 

Changing the notation to conform with the already estab- 
lished associated polynomial notation: 

An-i(Pi) = - (P3^+P2+* • ^ (l-27a) 

^n-2^Pi^ ^ ^PiP2‘*'PiP 3‘'’‘ • •■’■Pn-lPn^ ' Pi ° (l-27b) 

• • • 

Expressing equations 1-27 explicitly: 

Aj (Pi) = (-1) with: p^ = 0 and j = 1,2,..., n-1 

(1-28) 
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CHAPTER II 



Partial Fraction Expansion Using Associated Polynomials 
2.1: Introduction 

The basic associated polynomial theory developed in the 
preceding chapter has immediate application in the field of 
network analysis and network synthesis. Specifically, asso- 
ciated polynomials may be used as an alternate method for 
performing the Heaviside or partial fraction expansion of a 
ratio of polynomials. 

In network analysis, a method of finding the impulse re- 
sponse, h(t) , from the system's function is to perform a par- 
tial fraction expansion of the response's Laplace transform, 
H(s), and then use standard inverse transform tables. Many 
network synthesis procedures require that network functions 
be broken up into the sum of a niomber of terms which can be 
identified with network elements. In the case of the well 
known Cauer procedures, this involves partial fraction ex- 
pansions . 

The use of associated polynomials in partial fraction 
expansions is not the easiest method available when hand cal- 
culations are to be used. However, the calculations involved 
using associated polynomials are very adaptable to the dig- 
ital computer. Although an actual computer program is not 
presented here, this chapter presents the algorithm, or the 
step-by-step procedure for accomplishing a partial fraction 
expansion of any ratio of polynomials. 
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Consider the rational function: 



C(s) = 



m 

I b . 

i=0 ^ _ N (s) 

n . dTsT 

I 

i=0 ^ 



( 2 - 1 ) 



where a^^ = 1 and n ^ m. If equation 2-1 reoresents the trans- 
formed solution of an n order differential equation, then 
the denominator term is the characteristic equation of the 
system under study. It can be shown that the roots of the 
characteristic equation can occur in no more complicated form 
than a pair of complex conjugate roots or as multiple roots. 
Therefore, for convenience, partial fraction expansions using 
associated polynomials are developed using three different 
cases. These are: 

Case (1): C(s) contain first-order poles only. 

Case (2) : C(s) contains a pair of complex poles plus 

first-order poles. 

Case (3): C(s) contains multiple poles plus first order 

poles . 

2.2: Case (1) Partial Fraction Expansion 

If the roots of the denominator of equation 2-1 are each 
real and distinct, then by the fundamental theorem of algebra, 
the function, C(s), may be rewritten as: 



C(s) = 



N (s) 



(s-p^) ^s-p2) . . . (s-p^) 

The partial fraction expansion of equation 2-2 is: 

k 



( 2 - 2 ) 



^1 ^2 

C ( s ) = -j ^ + -j ^ + 

Ts Pj^ ) ( s P2 ) 



n 






(2-3) 
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The coefficient, k^, called the residue of C(s) at s = p^, 
i = l,2,...,n, is an undetermined coefficient of the expansion 
and is found by: 



(2-4) 

s=Pi 

Equation 2-4 is the classical method of finding the 
residues of a partial fraction expansion when all the roots 
are simple and distinct. This calculation is easy to per- 
form by hand calculations but requires several extra steps 
on a digital computer which will not divide by zero. Asso- 
ciated polynomial theory may be used instead of the above 
computations by the following method. 

Rearrange equation 2-4 as : 



k. = 
1 



N(s) 

DliT 



(s-p^) 



k. = 
1 



N(s) 

D(s)/(s-p^) 



i = 1 ,2 , . . . ,n 



(2-5) 



s=p. 



The denominator of equation 2-5 can be expressed using equa- 
tion 1-17 as: 



n 



= _l s^"^A. (p^) 



D(s) 

s-p. ~ "i 

3 = 1 



( 2 - 6 ) 



where A^ (p^) is the set of associated polynomials of D(s) 
evaluated at s = p^ . Hence equation 2-5 can be expressed as: 



k . = 

1 n 



N(P^) 



1=1 ^ 



i = 1 ,2 , . . . ,n (2-7) 



Example 2-1 



Expand C(s) = + 2 2s + 6 ^ partial fraction expansion. 

+ 3s^ + 2s 
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Solution : 



^,.1 _ N(s) _ 12s^+22s+6 _ ’^1 , ’^2 , ’'s 

- ora - - ^ 3 ^ 3 ^ 2 + 2 s ■ = =-^2 

Let: Pj^ = 0; p^ = -1; p^ = “2. 

Form the associated polynomials of the denominator, D(s): 

A3 (u) = 1 

A 2 (u) = 3 + u 
A j^(u) = 2 + 3u + 

Applying equation 2-7, the partial fraction coefficients are 

V = N(s= 0 ) _ 6 _ 

1 a]ToT+o 2 



where : 



and: 



Hence : 



where : 



and 



^2 = 



N(s= -1) 

A3^(-1) + (-1)A2(-1) + (-1)^A3(-1) 



Aj^(-l) = 0; A2(-1) = 2; A3 (-1) = 1 



N(s= -1) = 12 - 22 + 6 = -4 



k = ^ — 4 

^2 0 - 2+1 ^ 



^ N(s= -2) 

^ A^(-2) + (-2)A2(-2) + (-2)^A3(-2) 

Aj_(-2) = 0; A2(-2) = 1; A 3 (-2) = 1 

N(s = -2) = 48 - 44 + 6 = 10 
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Hence ; 



1 _ 10 _ r 

" 0+l(-2)+l(4) " ^ 

The partial fraction expansion of the function, C(s), is: 



rfs) - 12s +22S+6 _ 3 _ 4 5 

■ s'+3s^2s ■ = =-^1 * s+2 

2.3: Case (2) Partial Fraction Expansion 

When the denominator of equation 2-1 contains a pair of 

complex roots plus first order poles, then the function, C(s), 

is written: 



C(s) = 



N(s) 



(S-Pj^) (s-p^^*) (S-P2) . . . ^S-p^_^) 

The partial fraction expansion of equation 2-8 is: 



(2-8) 



Ic * 
^1 



n -1 



C ( S ) — / \ + ~7 apt" + ~i p + ... + I I — — ■■ 1 

(s-p^^ (s-pj^) (S-P 2 ) '®”Pn-l 



(2-9) 



As shown in reference 6 , the coefficients, and 
must be complex conjugates and either may be found by the 
same technique used in case ( 1 ) above. 

Example 2-2 

4 

Expand C(s) = 5 by a partial fraction expansion. 

s (s'^+ 2 s+ 2 ) 



Solution: 

C(s) 



N(s) ^ 4 = + 1^2 . 

s(s^+ 2 s+ 2 ) V® 



Let Pj^ = 0; P 2 = -1+j; ? 2 * = “l-j- 



Form the associated polynomials of the denominator, D(s): 

A^(u) = 1 

A2(u) = 2 + u 

A^ (u) = 2 + 2 u + u^ 
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Applying equation 2-7, the partial fraction coefficients are: 



N(0) _ 4 _ 
1 a^ToT+0 2 



2 



k 



2 “ 



N(-i+i) 

(-1+ j ) + (-1+ j) A 2 (-1+j ) + (-1+ j) ^A^ (-1+j ) 



Evaluating : 



Aj^(-l+j) = 2- 2 + 2j-2j = 0 
A2(-l+j) = 1 + j 



A3(-l+j) = 1 
N(-l+j) = 4 



hence : 



u -2 1 

^2 = T+J = 



j 



As the coefficient k 2 * must be the conjugate of k 2 , 



k * = 
^2 



-1-j 



The partial fraction expansion of the function, C(s), is: 
C (s) = — 



* = 2 + - 1-3 

s (s +20+2) 



= 2 ^ ^ ^ 

s S+I-3 s+l+j 



2.4: Case (3) Partial Fraction Expansion 

When the function, C(s) , is of the form: 



C(s) = 



N(s) 



K 



(s-pj^) (S-P2) . . . (s-p^_j^) 

The partial fraction expansion is: 



N(s) 
D (s) 



C(s) = 



m 



'"m-l 



, .K - .K-1 

(s-pj^) (s-pj^) 



s- 



Pl 



s-p. 



( 2 - 10 ) 



n-K 



s-p „ 
n-K 

( 2 - 11 ) 
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As shown by equation 2-11 the partial fraction expansion is 
modified to include as many additional coefficients as the 

order of the multiple root. 

(6 7 ) 

Several methods' ' are available for the evaluation 
of the coefficients of the multiple roots. Computer programs 
are also available ' ' that perform partial fraction expan- 
sions of expressions containing multiple roots but as shown 
in reference 7, there are several inherent sources of error 
that tend to degrade the numerical accuracy as the order of 
multiplicity of the root increases. 

This section introduces a procedure for expanding a mul- 
tiple root expansion, such as equation 2-10, using a combina- 
tion of associated polynomial theory and simple algebraic 
manipulations. The only source of error in the method is 
that the denominator polynomial, D(s), must be in factored 
form. 



The classical method 



( 10 ) 



of determining the partial frac- 



tion expansion of the function, C(s) is: 

N 



N . 
1 



C(s) = I I 



'iK 



i=l K=1 (s-p^) 



K 



( 2 - 12 ) 



where the coefficients, are given by: 



n = 1 d i 

iK (N.-K)l 

The interdependence of the successive derivatives in 
equation 2-13 is a possible cause of errors which cumulate 
as the multiplicity, N^, increases. A second method of 



-K r N. 1 

(s-p^) (s J 



(2-13) 



s=p . 
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using algebraic manipula- 



determining the coefficients, 
tions only is demonstrated in examples 2-3 and 2-4. 
Example 2-3 

Given the rational function : 



C(s) = 



ko 

2 Ti+TT 



(2-14) 



s(s+l)^ ° (s+1) 

Using the classical method of determining the coefficients, 
k^ and found to be: 



k^^ = s C(s) 



s=0 



k 2 = (s+1) C(s) 



= 1 



s=-l 



= 1/-1 = -1 



To find the coefficient, k^, multiply both sides of equation 
2-14 by the denominator of k^/ i.e., (s+1). Using the known 
values of k^^ and k 2 , combine the terms that do not go to zero 
as s-»— 1. There will be a cancellation of the undeterminate 
terms after the combination of terms. After the cancellation, 
let s = -1 and solve for k^* 



(s+1) 



s (s+1) 



= (s+l)-| + 



(s+1) k, 
(s+1) ' 



Let s-»— 1 and combine the non-zero terms: 

k. 



^^3 = 



1 

s ( s+1) 



Ti+^ 



(s+1) ko 
(s+1) 



but k 2 = -1, hence 

, _ (s+1) 

^3 - sTs+lT • 

Let s = -1, therefore the coefficient, k^, is: k^ = -1. 
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The denominator of the unknown coefficient will always 
be cancelled when the above procedure is used. The resulting 



partial fraction expansion is; 



C(s) = 



s (s+1) 



= ‘ (s4-i)2 ■ 



Example 2-4 



C(s) = 



s+2 



(s+1) (s+3) 



ki k^ k^ 



(2-15) 



As before: 



kj^ = (s+l)C(s) 



- 1+2 

s=-l (-1+3) 



1 

3 - ¥ 



and 



k^ = (s+3) C(s) 



s=-3 



-3+2 _ 1 
-3+1 2 



To 



find k^, multiply each side of equation 2-15 by (s+3) : 



(s+3) ^ (s+2) 
(s+3) ^ (s+1) 



(s+3)^k^ (s+3)^k2 (s+3)^k3 (s+3)^k^ 

(s+1) 3 2 (s+3) 



(s+3) 



(s+3) 



Combine all terms that do not go to zero as s->-3. 

k - s+2 _ Jl2 

^^3 " (s+3) (s+1) s+3 



but k 2 = hence 



, _ S+2-S/2-1/2 _ 1/2 (s+3) 

^3 “ (s+1) (s+3) " (s+1) (s+3) 



Therefore : 



^3 2(s+l) 



s=-3 



1 

■4 
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I 



The same procedure is used to find the coefficient, , i.e., 

(s+3) (s+2) (s+3)k^ (s+3) k2 (s+3) k^ (s+3) k^ 

(s+3)^(s+l) (s+3)^ (s+3)^ fs+3T 

Let s-^-3 and combine all non-zero terms. 



^4 = 



s+2 



(s+3) ^ (s+1) 



(s+3) 



But ^2 ~ ^ ^3 “ ~T' 



^4 = 



s+2-s/2-l/2+2s/4+3/4 _ 



s+3 



(s+3) 



(s+1) (s+3) 4(s+l)(s+3) 

Let s = -3, therefore the coefficient, k^ , is: 

The partial fraction expansion is: 
s+2 



k — — i. 
^4 8 



C(s) = 



- 1/8 1/2 

3 " (s+1) 



+ -1/4 ^ ;1Z 8 

(s+3)3 (s+ 3)2 Ti+3T 



(s+1) (s+3) 

In the preceding two examples, the algebraic manipula- 
tions required in the evaluation of the multiple root coeffi- 
cients are adaptable to the digital computer. The coefficients 
that were evaluated in the classical manner may be determined 
by associated polynomials, specifically by the use of Theorem 
1-2 and corollary 1-3. The combined operation is detailed 
below for two possible cases, i.e., when the repeated roots 
are of order two and three. 

Given the function: 
m 



C(s) = 



I b^s^ 
i=0 ^ 

(s-pj^) (S-P2) . . . 



(s-Pn-2> 



N(s) 

oTiT 



(2-16) 
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The partial fraction expansion of C(s) is: 

k . , k k, k „ 

C(s) + 1 + ... + T— 2^. 



(2-17) 



The coefficients, k^, k 2 , found using the 

procedure outlined in case (1) above. The multiple order 
coefficients, and k^ are found by the following method. 



m+1 



= (s-p ) C(s) 



m 



N(s) 



==Pm D(s)/(s-p^) 



(2-18) 



s=p 

^m 



The denominator of equation 2-18 can be expressed using equa- 
tion 1-19 as: 



D(s) h ? j-(h+2), , , 

~T = Z, s > 'A. (p ) 

(s-p )^ h=0 j=h+2 ^ 3 ^ 

*^m 

Evaluating at s = p^, the coefficient equals 

N(p„) 

, ni 



(2-19) 



m+1 n-2 , n • 

^ h ^ 3-(h+2)^ 

.hio j=h+2 



(2-20) 



In order to find the coefficient, k^, the algebraic method 
shown in examples 2-3 and 2-4 must be employed, that is, 
multiply each side of equation 2-17 by the denominator of k^. 



N 



%Wr ■ 



(s) 



m 



(S"Pm) 
^ in 



. m 

2 (s-p ) 

‘ m 



(s-p^) 



+ n-2 

• • • “T" / \ 

‘=-Pn-2'j 



(2-21a) 



Let s->p and combine each term that does not go to zero, 
'^m ^ 
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(2-21b) 



m+1 



N (s) [ (s-P^) <s-p2) . . . (s-p^_2) ] 

(s-Pni) (s-Pi) • • • (s-Pj^_2) 



s=p 



m 



The numerator of equation 2-21b will always contain a fac- 
torable term, (s-p ) . Since (s-p ) cancels in the combination 

'^m ^m 

of the two terms which individually become indeterminate as 
s->p , the combination is finite. Hence: 



N (s) (s“Pi) (S-P 2 ) ^®“Pn-2^ 



k = 
m 



(s-p ) 
m 



(S-pj^) (S-P 2 ) . . . (s-p^ t ) 



n-2' 



( 2 - 22 ) 



s=p. 



m 



Example 2-5 

Expand by partial fractions: 

C (~) - 1 I ^2 

■ - DTir - - * TiTIT ^ 



(s+1) 



Solution ; 



Let: p, =0; p^ = -1. 

1 •• m 



Form the associated polynomials of the denominator, D(s) 

(u) = 1 

A2(u) = 2 + u 
A^ (u) = 1 + 2u + u^ 

Applying equation 2-7, the coefficient is: 

1 Aj^(0)+0 ^ 

Applying equation 2-20, the coefficient, = k^ is: 

k = N(-l) 

3 •A3(-1) (-1) + [A2(-1) + (-1)A3(-1) ) 



^3 - =i+T-r ■ 
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Applying equation 2-22, the coefficient, = k 2 is: 

l-(-l)s 



^2 = 



s+1 



= -1 



s= -1 

The partial fraction expansion of C(s) is: 



C(s) = 



2 s 



1 . , 1 
2 ^ Ti+TT 



s(s+l)'^ ° (s+1) 

which agrees with example 2-3. 

When the expression, C(s), contains a third-order pole 
such that: 

m 

1 

C(s) = 5 = r z y (2-23) 

(s-p^) (s-pj^) (S-P2) . . . (s-p^_^) 

The expanded form is: 



r(~\ _ N(s) _ m+2 , ^m+1 , m , ^1 . . 

D(s) , v3 , ^2 (s-Pj^) (s-pj^) 



/ \ 3 , s2 

(s-p„) (s-p^) 



n-3 



(2-24) 

As above, the coefficients, k^, k 2 / • . . / found using 

the procedure of case (1). The multiple order coefficients, 
^m+2 ' ^m+1' ^m found by the following: 



'^m+2 



= (s-Pj,^) C(s) 



N(s) 



S=p„ D(s)/(s-p^) 



(2-25) 



®=Pm 



Using equation 1-21, the denominator of equation 2-25 can be 
expressed as : 

I [j-(h+l)]p 3-(h+2) ( j (2-26) 

(s-pj2 hil j=^2 ^ 3 - HI 
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Evaluating as s = the coefficient is: 

N (p_) 
m 



m+2 n-2 



I. p;*-V I 



(2-27) 



3 'm 



h=l j=h+2 

The coefficients, k ,, and k_ , are evaluated as in the 

m+1 in 

second-order example above. These coefficients are found to 
be : 

= <s-P2’ •• • *®‘Pn-3' 



(s-pj 



m+1 



(s-pj^) (S-P2) . . . 



S=D 



m 



(2-28) 



and 



N(s)-kj^^2<®-Pl> • <=-Pn-3>-’"m+l<^'Pm> •• • <®'Pn-3> 






(s-p ) 
- m 



(s-p^) (S-P 2 . . . (s-Pn-3^ 



S=D 



m 



(2-29) 



Example 2-6 

Expand by partial fractions: 



C(s) = 



s+2 



N(s) ^ 2i_ 
D (s ) s+1 



(s+3) 



(s+l)(s+3)-" """ (s+3)^ (s+3)^ 

Solution : 

Let: p^ = -1; p^ = -3 

Form the associated polynomials of the denominator, D(s) : 
D(s) = s"^ + lOs^ + 36s^ + 54s + 27 
A^(u) = 1 

A^(u) = 10 + u 
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A2(u) = 36 + lOu + 

A^(u) = 54 + 36u + lOu^ + 

Applying equation 2-7, the coefficient, k^, is; 

k = N(-l) 

^ Aj^(-1) + (-1)A2(-1) + (-1)^A3(-1) + (-1) 

where; N(-l) = 1; 

A^(-l) = 27; A2(-1) = 27; A3(-l) = 9; A^(-l) = 1. 



Hence ; 

k = 1 1 

^1 27-27+9-1 8 

Applying equation 2-27, the coefficient, k 2 , is; 

k = N(-3) 

2 (-3)A4(-3) + [A3 (-3)+2(-3)A^(-3) ] 

where ; 

N(-3) = -1; A3(-3) = 7; A^(-3) = 1. 



Hence ; 




-1 _ 1 
-3+7-6 2 



Applying equation 2-28, the coefficient, k 3 , is; 



k 



3 “ 



s+2-1/2 (s+1) 
s+3 
s+1 



1/2 



s+1 



s 



= -3 



k3 = -1/4 

Applying equation 2-29, the coefficient, k^ , is; 

s+2-1/2 (s+D+1/4 (s+1) (s+3) 

_ (s+3)^ ^ 1 

4 s+1 ~ 4 (s+1) 

k4 = -1/8 



-3 
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CHAPTER III 



Polynomial Sensitivity Theory 
3.1: Introduction 

In this chapter several simple algebraic operations in- 
volving polynomial sensitivity are investigated. Specif- 
ically, associated polynomial theory is used to study the 
relationships between variations of polynomial coefficients 
and polynomial roots. Both large and incremental variations 
are considered. As is shown in the next two chapters, the 
results derived from this section have very useful applica- 
tions, particularly in the field of sensitivity analysis 
(Chapter Four) and in self-adaptive control systems (Chapter 
Five). In order to obtain these results, a basic ortho- 
normality relationship between the set of associated polyno- 
mials and the roots of the defining polynomial is defined. 

In addition to being used in the derivation of sensitivity 
expressions, the orthonormality relationship can be used to 
solve a variety. of linear problems that involve matrix in- 
version. 

3.2: Coefficient Sensitivity 

The effect of variations of a root on the coefficients 
of a polynomial can be expressed by the following theorem: 
THEOREM 3-1 

The variations of the coefficients of a polynomial with 
respect to the variation of a root may be determined using 
polynomials by the relationship: 



This result holds for large variations when is a simple 
root. 

Proof : 

Let a monic polynomial, F(s) , of degree n with n distinct 
roots be expressed as: 



n 



n . _ „n-l ^ , t 1 1 

s + . . . + a^^ = I a^S"^ ; a^ = 1 



F(s) = s“ + a ,s 

" ^ “ jio ' 

(3-1) 

By the fundamental theorem of Algebra, the polynomial may 
also be written as: 
n 






w • • /Pv^) 



n 



(3-2) 



Let Ap^ denote a change in the root, p^^. The polynomial, 
f (p, ,p^ , . . . ,p. +Ap . , . . . ,p ), may be expressed in a Taylor se- 
ries expansion about the point (p^^ ,P 2 , . . . ,p^^ , . . . ,p^^) by: 

f (p^,P2,. . .Pj^+Ap^,. . .p^) = f (p^,p2, . . . /P^/ . . . /Pj^) + 

3f(P'j^,P2/' •• /Pj^/..* 'Pj^^ 



3P) 



APj_ + 



3 f(Pi /Pn/*** /P: /Pn^ 9 

1/2 ^ ^ 2- (Ap.)^ + 



3p. 



1 3^f(p^,p2,... ,p^,...,Pj^) 

HI r-H <^Pi^ 

3p. 



(3-3) 



However, from the factored form given in equation 3-2: 
3f(pj^,P2/* /?£/• •• / Pj^^ 



3p. 



F(s) 

(s-p^) 



(3-4) 
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and; 



9 



9P, 



n 



= 0 : n I 



(3-5) 



Equation 3-5 applies when p^ is a simple root only. Hence 



equation 3-3 becomes 



F(s) 



f (p^,P2,. . . ,Pj_+Ap^,. .. ,Pj^) = f(pj^,P2/. .. ,Pj^) - Ap^ 



Pi 
(3-6) 



and ; 



AF = -Ap. 



F(s) 
1 s-Pi 



(3-7) 



where ; 

AF = f(Pj^fP2/***^Pj_"^Ap£/««.fPj^) — f(p2^/P2/«»«/ Pj^ ) ( 3 — 8 ) 

The expression AF in equation 3-8 may be expressed in terms 
of the changes in the coefficient, Aaj^, of the original 
polynomial, F(s) . Note that: 

f (pj^,P2, . . . ,Pj^+Ap^, . . . ,Pj^) = (s-pj^) — (s-p^+Ap^) ... (s-p^ 



(3-9) 

Sxibstituting equations 3-2 and 3-9 into equation 3-8 and s\ib 
tracting the coefficients of equal power of s reveals that; 

n 

Aa. 

j=l 

Using Theorem 1-2, i.e.,; 



= 'i Aa. 

T-1 



(3-10) 



F(s) 

S-p. 



I A. (p.)s^"^ 
j=l 3 1 



(1-17) 
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Equation 3-7 is rewritten as: 



AF = -Ap. y A. (p. ) 
^1 11 
3=1 -> 

Equating equations 3-10 and 3-11: 



(3-11) 




(3-12) 



Equating the coefficients of equation 3-12, it follows that: 



equation 3-13 holds in the large, i.e., for arbitrary changes 
in Ap^, when the root is simple. Similar precise equations 
may be derived for repeated roots. Equation 3-13 is useful 
for determining the change in the coefficients of a polyno- 
mial when there is a change in a simple root. Equation 3-13 
cannot be used to determine the change in root locations 
caused by changes in any coefficient of a polynomial since 
a change in general can be expected to change the location 
of several roots. Thus the inverse of equation 3-13 is not 
meaningful . 

Corollary 3-1 

Taking the limit of equation 3-13 as the change in the 
root, p^ , goes to zero, then: 





= -Aj (p^) : i, j = 1,2 



(3-14) 
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(s+1) (s+2) (s+3) 



Example 3-1 
Given the polynomial : 

F(s) = + 6s^ + 11s + 6 = 

The associated polynomials are: 

A^Cu) = 1 

A 2 (u) = 6 + u 

Aj^(u) = 11 + 6u + 

Let the roots of F(s) be defined as: 

Pi = -1; P 2 = -2; P 2 = -3 

Assume p^^ is decreased by 8 to a new value p^ + Ap^ = -9. 

Hence : 

Ap^ = -8 

The original polynomial is changed to: 

Fj^(s) = s^ + 14s^ + 51s + 54 = (s+9) (s+2) (s+3) 

By inspection, the change in the coefficients are seen to be: 
ASq = 48; Aa^ = 40; Aa2 = 8; Aa^ = 0 

In order to verify theorem 3-1 numerically, the changes 
in the coefficients can also be found by use of equation 3-13. 
For i = 1 : 

Let j = 3, or: 

Aa2 = -A2(pj^)Apj^ 

numerically : 

Aa2 = -d) (-8) = 8; 

which agrees. 

Let j = 2, or: 

Aaj^ = -A 2 (pj^) Ap^ 
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numerically : 



Aa^ = -[6+(-D ] (-8) = 40; 



which agrees. 



Let j = 1, or; 



= -A^(pj^)Apj^ 



numerically : 



AaQ = -[11+6(-1) + (-1)2] (-8) = 48; 



which agrees. 

3.3; The Basic Orthonormality Relation 

In order to derive the inverse of Theorem 3-1, that is, 
the variation of a root of a polynomial with the variation 
of its coefficients, the orthonormality relationship of asso- 
ciated polynomials must first be derived. Prior to this 
derivation. Theorem 3-2 is introduced. 

THEOREM 3-2 

The derivative of a polynomial may also be found from 
its associated polynomials by; 




I u^“^A. (u) 
j=l ^ 



n 



Proof; 



It was shown in the proof of Theorem 1-2 that; 



F(s) -F(u ) 
s-u 




(1-14) 



Let s->-u, then 




II * ^ 

= lim ^ s^ "^A . (u) 

-1 = 1 3 



s->-u j=l 



(3-15) 
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or 



F' (u) = J . (u) 

j = l ^ 



(3-16) 



Example 3-2 

Given the polynomial : 

F(s) = + lOs^ + 35s^ + 50s + 24 = (s+1) (s+2) (s+3) (s+4) 

t 

The associated polynomials are: 

A^(u) = 1 

A^Cu) = 10 + u 

A2(u) = 35 + lOu + 

Aj^(u) = 50 + 35u + lOu^ + 

The derivative of the polynomial, F(s), is expressable in 
two ways: 

(1) : Classical Method: 



= 4s^ + 30s^ + 70s + 50 
ds 

(2) : Associated polynomial method; 

Applying Theorem 3-2; 

= [Aj^(u)+A2(u)u+A2(u)u^+A^(u)u^]^_g 

= 50 + 35u + lOu^ + 

+ 35u + lOu^ + 

2 3 

+ lOu + u-^ 



dF(s) 

ds 



50 + 70u + 



30u^ + 



4u' 



u=s 
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THEOREM 3-3 



The basic orthonormality relationship between the set of 
associated polynomials and the roots of the defining polyno- 
mial is : 



where is the Kronecker delta function which vanishes ex- 

cept when the indices are equal. 

Proof : 

The derivative of a polynomial is expressable using 
Theorem 3-2 as : 



then it can be shown that, if p^ is a simple root of F(s) : 




, . . . ,n 



F ' (u) = ^ u^ . (u) 

j=l 3 



(3-16) 





(3-17) 



From equations 3-16 and 3-17, it follows that: 




when i ^ k 
when i = k 



(3-18) 



Equation 3-18 may be rewritten: 




ik 



i,k = 1,2 



, . . . , n 



(3-19) 



where 6., is the Kronecker delta function. 
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(3-20a) 



Let [a] and [6] denote the matrices who elements are: 
A. (P,) 



a . . = 



- 1 1 



and 



ij f' (Pj^) 
j -1 






^ — if2/««.;nj k — 1^2f«««/n (3 20b) 

If I denotes the identity matrix, then equation 3-19 may 
be written: 

[a] [6] = ^ij^jk = ‘^ik = ^ 

Since is a Vandermonde matrix, it is non-singular and: 



n 



IB! [a! = 6 ,^. c. . = I 



(3-22) 



Hence : 



k -1 



n A^(Pj)Pi _ 
F' (P^) 



(3-23) 



^,k l, 2 ,»**,n 

Theorem 3-3 may also be verified by the following example: 
Example 3-3 

F(s) = s^ + 6 s^ + 11s + 6 = (s+1) (s+2) (s+3) 



Let: 



Pi = -1; P 2 = -2; P 3 = -3. 

The associated polynomials of the given polynomial are: 

A3(u) = 1 

A2(u) = 6 + u 
A^(u) = 11 + 6 u + u^ 

The derivatives of the polynomial evaluated at the three 
roots are found using equation 3-16 to be: 
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F'(p.) = A^(p.) + A 2 (p.)p. + A3(P.)p.2 
F' (p^) = 6 + 5(-l) +1=2 
F' (P 2 ) =3 + 4 (-2) + 4 = -1 
F' (P 3 ) =2 + 3 (-3) +9 = 2 



Forming the [a] and [3] matrices according to equation 3-20 



Aj^(pj^) 


A2(Pi) 


A3(pi) 




F' (p^) 


F' (pj^) 


F' (p^) 






Ai(p2) 


A2 (P 2 ) 






6/2 5/2 1/2 

3/-1 4/-1 1/-1 


f'TpJT 


FlpP" 


pn:pjT 




Ai(p3) 


A2 (P3) 


A2(P3> 




2/2 3/2 1/2 




FTp?" 









3 


5/2 


1/2 




-3 


-4 


-1 




_1 


3/2 


1 / 2 . 





0 

Pi 


0 

P 2 


0 

P 3 




1 — ! 

1 


1 


1 




1 

Pi 


1 

P 2 


1 

P 3 


= 


-1 


-2 


-3 




2 

[Pi 


2 

P 2 


2 

P3_ 




_^1 


4 





It can be shown that: 

[a. .] [3 ■!,] = I (3-21) 

Likewise : 

[&.,] [a. .] = I (3-22) 

J ^ J 
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Expanding the terms in equation 3-22; 



0 


0 


0 


Pi 


P2 


P3 


1 


1 


1 


Pi 


P2 


P3 


2 


2 


2 


LPi 


P2 


P3J 



(p^) 


^2 ^Pl^ 


A 3 (Pi) 


F' (p^) 


F' (p^) 


F' (Pj^) 


Ai(p2 ) 


^2 ^P2^ 


A3 (P2) 


FTip- 


F-(p2) 


FTi?" 


Ai(P 3 > 


^2 ^P 3 ^ 


A3 (P3) 


_fTp 7 




F'(p 3 >_ 



= I 



Hence ; 



0 . _ 0 , _ 0 
Pi F' (p^) P2 F^TpJT P 3 F* (P 3 ) 



= 1 



(3- 



Or in matrix notation, 



^^j=l,k^ ^“i, j=l^ “ ^ 



(3- 



likewise : 



[6 



j=2,k^ ^“i,j=2^ “ ^ 



(3- 



and; 



^^j=3,k^ ^“i, j^3^ " ^ 



(3- 



In addition, the following relationships also hold: 



I3j^l,k] 


= 0 


(3 


f^j=l,k^ ^“i,j = 3^ 
« 

• 


= 0 


(3 


« 

>^3=3,k''“io-2l 


= 0 


(3 



The above matrix products, equations 3-24 and 3-25 



be expressed explicitly as stated in Theorem 3-3, i.e.. 



i=l ^Pi^ '^’^3 



(3 



24a) 

24b) 

24c) 

24d) 

25a) 

25b) 

25c) 

can 

23) 
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The orthonormality relation given by equation 3-23 is a 
basic result which permits one to solve a variety of linear 
problems. As has been shown in references 1 and 11, this 
orthonormality relation permits the theoretical investiga- 
tion of inverse matrix problems, such as finding the inverse 
to the Vandermonde Matrix in a form suited for numerical cal- 
culations . 

3.4: Root Sensitivity 

The effect of variations of the coefficients of a polyno- 
mial on any particular root is expressable by the following 
theorem : 



THEOREM 3-4 



The variation of the roots of a polynomial with respect 
to the variations of a coefficient of the defining polyno- 
mial is determined by the relationship: 



8 5i • 

3 



Pj^ ] — 0,l,...,nl 

’FTp.) ‘ , 



Proof : 

Corollary 3-1 established that: 
9a . 



-A. (p.) : i,j = l,2,...,n 

dp^ 3 1 

In addition the orthonormality relation of Theorem 3-3 is: 



n , / V k-1 



= 5 



i=l 



F' (p^) 



kj 



(3-23) 



Let [y] and [n] denote the matrices whose elements are: 

— Aj(p^) r — 1,2,. ..,n^ j — 1,2,. ..,n (3— 26a) 
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and 



k-1 



1 ^ f 2 f » 9 » f ^ f 2 f » • • f n 



(3-26b) 



'^ik F' (p^) 

Using the matrix notation introduced above, equation 3-23 



can be rewritten: 



[y] [n] = I 



(3-27) 



Hence, the matrix, y, (whose elements are A.(p^)) is the in- 
k-1 

verse of the matrix, r\ , (whose elements are /F'(p^)) 

given i,j,k = l,2,...,n. In a similar manner, it has been 

shown that the matrix with elements 9p./9a. , is the in- 

j-l 

verse of the matrix whose elements are 9aj_j^/9p^, that is, 
if : 



and 



[A] = 



[B] = 



9a.-i 



9a 



k-1 



9p, 



Then: 

[A] [B] = I 

but from corollary 3-1: 

[B] = -[y] 

Hence, from equations 3-27 and 3-28d: 

[A] = -[n] 



(3-28a) 



(3-28b) 



(3-28c) 



(3-28d) 



(3-28e) 



Writing out the terms of the matrices, as given by equations 
3-28a and 3-26b yields: 
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9P]_ 


9pj_ 


9Pl 




1 

O 


1 

Pi 


n-1 

Pi 


9a 

o 

9p2 


9a^_ "• 


• 




F' (p^) 

0 

P2 


F' (P3_) 


•*'F'(P3_) 


9ao 

• 




• 


— — 


• 




• 


• 

• 

^Pn 




^Pn 




• 

• 

0 

Pn 




n-1 

Pn 




• • • 








• • • 





Thus in general : 



9Pi i = 1,2 , . . . ,n 

9a^ F (p^) j ^ 0,l,...,n-l 

Corollary 3-2 

The total change in the value of any root of an n 
order polynomial due to an incremental change in any or all 
coefficients of the polynomial is found by: 



(3-30) 



-th 



dp. = 



9p. 



9p. 



^ da,, + — — da, + 

9aQ 0 9a^^ 1 



9p. 



9a 



da 



n-1 



h-1 



where 



9p. ^ 
9a 

m 

Substituting and factoring: 

‘o 



m 

Pi 

F' (p^) 



(3-30) 



dp. = 



da„ + p.da, + 
n 1 



+ p . da 

•^1 



n-1 



-F' (p^) 



(3-31) 



3.5: Root Solving Using Associated Polynomials 

Equation 3-31 of Corollary 3-2 may be used in a root 
solving process if a polynomial with known simple roots can 
be found near the unknown polynomial. Since equation 3-31 
deals with differentials, the process is limited to 
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polynomials whose roots are close in value to the known polyno 

mial, i.e., small variations in the coefficients between the 

1 

two. The process is summarized as; 

(1) : Find a polynomial of equal order with known, non- 

repeated, roots, such that the difference in the coefficient 
values are "small". The accuracy of this root solving proc- 
ess is determined by the value of the variations in the coeffi 
cients . 

(2) : Using the roots of the known polynomial and the 

values of the change between the coefficients of the known 
and unknown polynomials, apply equation 3-31. 

(3) : The change in the known roots determines the 

values of the roots of the unknown polynomial. 

The following examples demonstrate how this root solv- 
ing process may be employed. 

Example 3-4 

The roots of a quadratic polynomial are easily deter- 
mined by the quadratic formula. However, such analytical 
expressions, in general, do not exist for higher order poly- 
nomials. This example will use a quadratic polynomial as 
an illustration of the root solving method available with 
associated polynomials. 

Given the quadratic polynomial: 

G(s) = s^ + 4.7s + 5.3 

In order to find the roots of G(s) , note that the polyno- 
mial , F (s) , is : 

F(s) = s^ + 5s + 6 = (s+2) (s+3) 
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The variations in the coefficients of F(s) from the given 
polynomial, G(s), are small enough such that the process may 
be employed with negligible errors. 

Let the roots of F(s) be defined as: 

“ “2; P 2 — “3 . 

Evaluating the denominator of equation 3-31: 

F' (s) = 2s + 5 



Therefore : 

-F' (pj^) = -1; -F' (P 2 ) = +1 

The variations of the coefficients are identified as: 
daQ = -.7; da^ = -.3; da2 = 0. 



The change in the root, p^^, of F(s) 

-.7+(-2) (-.3) 



Hence ; 



api = ^ 

Pj^ of G(s) is -2 + .1 = -1.9 



The change in the root, P 2 , of P(s) 

-.7+(-3) (-.3) 



dp, = 



+ 1 



is : 

= .1 



is : 



= .2 



Hence: P 2 of G(s) is -3 + .2 = -2.8. 

The factors of the polynomial, G(s) , are: 

(s+1.9) (s+2.8) 

As a check on the above results : 

(s+1.9) (s+2.8) = s^ + 4.7s + 5.32 

Example 3-5 

Given the third order polynomial: 

G(s) = s^ + 5.85s^ + 10.44s + 5.39 
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To find the roots of G(s), note that the polynomial, F(s), 



IS : 



F(s) = s^ + 6s^ + 11s + 6 = (s+1) (s+2) (s+3) 



Let the roots of F(s) be defined as; 

~ “1? ?2 “ “2; = “3. 

The denominator of equation 3-31 is; 



Therefore 



F' (s) = 3s^ + 12s + 11 

F' (pj^) = 3 - 12 + 11 = 2 
F' (P 2 ) = 12 - 24 + 11 = -1 
F' (P 3 ) = 27 - 36 + 11 = +2 



The variations of the coefficients, [a^ of G(s) - a^ of F(s)], 
are ; 

daQ = -.61; da^ = -.56; da2 = -.15 
The change in the root, p^, of F(s) is; 



daj+pj^daj^+pj^ daj 




^ (-.6D-K-1) (-.56) + (-l)^(-.15 ) ^ , 

api _2 .X 



Hence 



; pj^ of G(s) is -1 + .1 = -.9. 



The change in the root, P2 , of F(s) is 



dai+P2dai+p2 da2 

'*P2 = -F' (p,) 



dp. 



= (-.61) + (-2) (-.56)+ (-2)^ (-.15) ^ 



+ 1 



09 
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-2.09. 



Hence: P 2 of G(s) is -2 - .09 = 

The change in the root, p^, of F' (s) is: 

dao+P3dai+p32da2 

^P3 = -F'(p3) 

dp^ = (--eP + C-S) (-^56) + (-3)^(-.15) ^ 

Hence: P 3 of G(s) is -3 + .14 = -2.86. 

The factors of the polynomial, G(s), are: 

G(s) = (s+. 9) (s+2 . 09) (s+2 . 86 ) 

As a check on the above results : 

(s+. 9) (s+2. 09) (s+2. 86 ) = s^ + 5.85s^ + 10.43s + 5.38 
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CHAPTER IV 



Sensitivity Coefficients 
4.1: Introduction 

In the analysis of a system whose motion can be described 
by a set of differential equations',- in addition to obtaining 
solutions as functions of independent variables, the engineer 
would like a knowledge of the variations of the solutions 
with respect to the parameters of the system. Since it is 
impossible to realize a physical system that is identical to 
a mathematical model, the influence which parameter varia- 
tions in the physical system have on the system's behavior 
is an important condition for reliable design. 

Variations in a circuit element will cause a disolace- 

ment of the poles and zeros of driving-point and transfer 

functions. As a direct consequence, these variations and 

resulting displacements Effect the network response. One 

method of approaching this problem is to plot the displaced 

poles and zeros and then determine the new response. As this 

( 12 ) 

method involves considerable work, it has been shown that 
the work involved may be reduced by testing the driving-point 
or transfer function to determine incremental changes in a 
particular parameter. This method is called "Sensitivity 
Analysis" . 

(13) 

Tomovic introduced sensitivity coefficients for 
error and sensitivity analysis of dynamic systems. He has 
shown that determination of the influence of parameter tol- 
erances on the behavior is based on a knowledge of these 
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coefficients. They show the effects of incremental parameter 
variations on the system's behavior by giving the variation 
of the system response for any given input. 

Systems which contain parameters that undergo large 
variations are studied in the next chapter. This chapter ex- 
presses Tomovic' sensitivity coefficients analytically in 
terms of associated polynomials. Tomovic' work demonstrated 
how these coefficients may be obtained from a computer model 
and analytic expressions for them have not been available. 
Applications of sensitivity coefficients for error and stabil- 
ity analysis has been amply covered in reference 13 and is 
not covered here. 

4-2: Theory of Sensitivity Coefficients 

Any dynamic system may be represented by an input, an 
output and a transfer function relating the two as shown in 
figure 4-1. 



Input y(t) 


System 


Output x(t) 









Figure 4-1; Dynamic System Basic Block Diagram 

If the Laplace transformation of the input signal, Y(s), 
is equal to one, then the response of the system is dependent 
on n parameters and is defined by; 



b,.+b, s+. . .+b s 
u 1 n 



n 



X(s,aQ,aj^ 



/ • • • / 




n 



(4-1) 




+anS 
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Tomovic introduced the complex sensitivity coefficient, 



V^{s), which can be defined with respect to the relative 



variation of the parameter a^ , as: 

3X(s,aQ,a2^,» • • , 



V^(s) = 



8 In a^^ 



(4-2) 



It follows from the Laplace integral's property of 
differentiability with respect to the parameters that the 
real sensitivity coefficient, Vj^(t), can be obtained by the 
inverse transformation: 



v^(t) 



3X(t,aQ,aj^,» ♦ ♦ / 
9 In a . 



(4-3) 



1 

In his text, Tomovic continues from this definition to 
a method of determining the sensitivity coefficients by 
modeling on an analog computer. This section uses asso- 
ciated polynomial theory for the determination of the sensi- 
tivity coefficients in terms of the coefficients of the 
system's transfer function which therefore makes them amen- 
able to digital computer computation. 

The Laplace transformation of the response, equation 4-1 
may be written as : 






n . n . 

1 b.s^ ,_1. \ b.s^ 

Jo 

where the denominator has been factored. Expanding equation 
4-4 into its partial fractions: 

k, k., k 



X ( S , a • ) = i' I + y r 

l' (S + Pj_) (S+P2) 



+ • • • + 



n 



(4-5) 
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From equation 4-5, two operations may be noted; 

( 1 ) : 



9X(s,a^) 



k. 

1 



( 2 ) ; 



9p, 



9X(s,a^) 

~W. 

1 



(s+p^) 



Ti+pT]" 



(4-6) 



(4-7) 



The response is a function of the parameters a^ and the 
residues , k . . 

The complex sensitivity coefficient, equation 4-2, may 
be written as ; 

9X(s,a. ) 






n 



= a. I 

^ jil 



ax(s,a^) 9pj 9X(s,a^) 9k^ 

9p. aa."*” 9k. 9a. 

^3 X 3 X 



(4-8a) 



(4-8b) 



From associated polynomial theory, the variation of a 
particular root with the variation of a coefficient is ex- 
pressable as: 



9 p^ Pj^ i — l,2,..«,n 

9a . ~ F' (p. ) . r\ 1 

] 3 = 0 ,l,...,n 



(4-9) 



For notational purposes, let: 



9aT 

1 



9P2 9Pn 

^1' 9a. ^2' •**' 9a. "n 

1 1 



(4-10) 



Using equations 4-10, 4-6 and 4-7, the complex sensi- 
tivity coefficients, equation 4-8 may be simplified to: 



V^(s) = I 

j=l 



-k. - 9k. 

3 0 + — 1 

/ . \ 2 ^ j (s+p.) 9a. 

(s+Pj) 3 1 



(4-11) 
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Equation 4-11 may be further simplified by noting: 
3k. 3k. 3p. 



3a. 3p . 3a . 
1 1 



/ J 1,2;. ../n 



(4-12a) 



or using equation 4-10/ 
3k. 3k. 



3T- ' a?' ' 3 ' " 

1 ^3 



(4-12b) 



The residue, k^ , in equation 4-12b may be expressed in terms 
of the original equation for the Laplace transformation of 
the response. From that equation, for notational purposes 
let : 



n 



N (s) = ^ b 



3 



j = 0 



n 



D(s) = J a.s^ 
j = 0 => 



(4-13a) 



(4-13b) 



Form the set of associated polynomials of the denominator 
D(s). From associated polynomial theory, the residue kj is 
expressable as; 



k . = — 
3 n 



N (p.^) 



y (p.)’^ (p.) 

mil 3 3 



3 ~ 1/2,. ..,n 



so that; 



3k . 



i = ^ 






9p . dp n T 

^ ^ I (p,)"'‘V(p.) 

' mil 3 "■ 3 



J 



(4-14) 



(4-15) 



Using equations 4-12 and 4-15, the complex sensitivity coeffi- 
cient equation 4-11 may be written as: 
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V^(s) 



Or 



V^(s) 



a . 

1 



®1 / ’^1 ^ 3 



N(Pi) 



TiT5i) s+p 3p^; n , 

I, Pi %<Pl>. 

m=l 



+ 



^2 I ^2 



N(Po) N 



Ts+pTrl s+p, n 



m=l 



/ 



= a . 



=n ( 


-k 

" + 


9 


N(p ) ^ 


*=-*-Pn* \ 


s+Pn 


9p ‘ 
^n 


' n , 

'I p A (p ) 

m ny 


5 Qj 


-k . 


^ 9 


N(p.) 


ii^s+p.) 


s+p . 

L : 


9P^ ^ 1 

1 (p.) 

m=l ^ 3 . 



(4-16a) 



(4-16b) 



This expression is cumbersome but does give an analytic ex- 
pression for the transform of a sensitivity function. 
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CHAPTER V 



A Two Parameter Adaptive System Using Associated Polynomials 
5.1: Introduction 

A feedback control system can be represented by the basic 
diagram as shown in figure 5-1, or by some modification or 
extension to it. The plant, control and feedback systems are 
composed of real physical components . For purposes of anal- 
ysis and design the physical system is usually modeled math- 
ematically in order that the Laplace Transforms may be em- 
ployed. 




Figure 5-1. Feedback Control System 

The characteristic equation of the system is expressable as: 

1 + HG(s) = 0 

where : 

H(s) is the feedback transfer function 
G(s) is the direct or forward transfer function com- 
posed of the individual transfer functions of both the control 
system and the plant. 

HG(s) is the open-loop transfer function. 

The characteristic equation determines the general form 
of the system's natural response and its roots give the time 
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constants of the transient exponential decay factors. Stab- 
ility of the overall system is found by factoring the char- 
acteristic equation. The characteristic equation has real 
coefficients which are algebraic combinations of physical 
parameters , which are real numbers . 

Basic feedback control theory usually assumes time- 
invarient elements for purposes of design and analysis. How- 
ever parameters of the system may be subject to variations 
due to a changing environment or similar phenomena. If there 
is a constant loop gain or a fixed compensation scheme and 
a system parameter were to vary, the closed looo performance 
could be seriously deteriorated and even system instability 
occur. Examples of such a changing parameter may be found 
in different operations such as the change in the aero- 
dynamic parameters of an aircraft as it changes altitudes 
or the effect of gravity on system parameters in the Ship- 
board Inertial Navigation System, (SINS) , of a deep diving 
submarine . 

If parameter variations are large and slow compared 
with the normal system's response times, it is possible to 
design for the capability of continuously measuring the 
variations and then changing compensation so that the sys- 
tem's performance criteria are always satisfied. Such sys- 
tems which automatically exhibit the characteristics of 
adapting to variations in system parameters, are called 
Self-Adaptive control systems. 
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Figure 5-2. Control System with Variable Parameter 6 

The parameter B of the plant shown in figure 5-2, is 
known to vary slowly with time or environment.* This param- 
eter appears algebraically in the coefficients of the sys- 
tem's characteristic equation. Since the dominant root val- 
ues are determined by the coefficients of the defining equa- 
tion, the system's dominant pole location is a function of 
the value of the variable parameter B. As the value of B 
varies so also will the location of the dominant root or 
roots of the characteristic equation hence the system's re- 
sponse to an input and decay factors will change. If the 
new value of B the location of the dominant root or 

roots into the right-hand s-plane, the system becomes un- 
stable . 

A block diagram of an adaptive control system is shown 
in figure 5-3. As before the plant has a parameter B which 
is known to vary with time or environment. The "identifica- 
tion and adapter" block continuously measures the input m(t) , 

*In the analysis which follows, the assumption of quasi- 
stationary responses is made in which parameters are assumed 
to be varying slowly such that the response is approximately 
given by the system's response with the fixed parameters re- 
placed by a time varying function. (23) 
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Figure 5-3. Self-Adaptive Control System 

and the output of the plant in order to identify the param- 
eter B. When a change in 3 has been detected the "identifica- 
tion and adapter" block adjusts the variable compensation 
block parameter a such that the system is returned to its 
original specification. The design of the "identification 
and adapter" block is the critical portion of any self- 
adaptive system and involves two requirements: 

(1) Some method must be devised for identifying the 
change within the system's plant, 

(2) A means of adjustment must be provided which is 
capable of counteracting the effects of the parameters varia- 
tions . 

Identification of any change in the parameters of a 
plant have been covered in references 14 and 15. References 
16 and 17 are concerned with methods for providing adjustment 
schemes. This chapter shows how associated polynomials 
theory may be used in the development of a "compensating 
ratio" between two variables parameters of a feedback con- 
trol system. It shows that a change in a variable parameter 
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in the plant may be compensated for by a change in another 
variable parameter such that dominant root performance or be- 
havior is controlled, 

5.2: Basic Theory 

The transfer function of a simple feedback control sys- 
tem as shown in figure 5-1 is written as: 

■"'=> <5-i) 

The denominator of equation 5-1 when equated to zero is the 
characteristic equation of the system. This equation is a 
polynomial of degree n and may be represented as: 

W(s) = a s^ + a - s^ ^ + ... + a^ = 0 (5-2) 

n n-1 0 

By the fundamental theorem of algebra, the characteristic 
equation has n roots. Hence, equation 5-2 may be written: 
W(s) = (s+p^) (S+P 2 ) . . . (s+Pj^) = 0 (5-3) 

Consider the coefficients of the characteristic equa- 
tion, i.e. , a^, a^, ..., a^. Assume any two variable param- 
eters a and B appear linearly in any number of coefficients, 
i.e.. 



a. = b.a + c-B + d. 
D D D D 



(5-4) 



: ^ 0,1,..., n 



where b., c., and d. are real constants. Note that the 
3 3 3 

coefficient a^ is a linear function of a and B. This repre- 
sentation of variable parameters in the coefficients of the 

characteristic equation was first introduced by D. D. Siljak 

(19) 

and later expanded by G. J. Thaler, et al. , m the 
development of parameter plane theory. 



( 18 ) 



66 



From general polynomial theory, the change in any root is 
determined by the change in the coefficients of the defining 
polynomial and vice-versa. The roots of the characteristic 
equation can be considered as functions of the characteristic 
equation's coefficients. Hence the differential p^ (aQ ,a^ , . . . ,a^^) 
can be expressed as : 



9p. 9p. 

dp^ = ^ ^ da^ + — ^a, + 



9ao 0 



9^ '‘I 



9Pi 

+ ^-ida 
9a n 
n 



(5-5) 



From equation 5-4, the coefficients, a^: j = 0,1,. ..,n, 
are themselves functions of the variable parameters a and 3. 
Hence the differential of any coefficient is; 



9a. 9a. 

Siabstituting equation 5-6 into equation 5-5: 

9a, 



(5-6) 



3p. r ^ 3p. 



9 3.^ 



9a 



-da + 



93 



^3 



9a 



9p. 9 a 

+ ^ 1.^-^a + 777^33 



9a 9a 
n 



93 



(5-7) 



From equation 5-4 : 

9a . 



v- ^ = b . 
9a j 



9a 



93 



(j — 0,1,. ..,n) 



(5-8) 



= c 



Substituting equation 5-8 into equation 5-7 and then factor- 
ing da and d3; 
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( 5 - 9 ) 



r 3p^ 9p^ 3 pri 

1 L 0 9aQ 1 9a^ n 9aJ 

r ^Pi ^Pi ^Pil 

9a^ ^1 9aj_ +•••+ ‘^n 9a^J 



dB 



‘0 " "“1 

From associated polynomial theory, the variation of a polyno- 
mial's roots with the variation of its coefficients is exores- 
sable as : 

j 



!£i 

9a. 



Pi 
F' (p. ) 



j — 0,1,. ..,n 



(5-10) 



Using the above equation, the differential of any root, , 
may be written as ; 



n 



b +b p +...+bnPi 
dp^ = sT-TT^ — ^ da + 



-F' (P^) 



Crt+C, p . .+c p 
0 1^1 n'^i 



n 



-F' (p.) 



dB 



(5-11) 



In this chapter the term "compensating ratio" is used 
for the mathematical relationship between two variable param- 
eters such that the system's dominant performance or behavior 
is controllable. Since the dominant root of a feedback con- 
trol system may be either real or a complex conjugate pair, 
four different cases of compensating adaption are possible. 
These cases are: 

(1) The dominant root is real and is to be held con- 
stant, 

(2) The dominant roots are a complex conjugate pair 
and are to be held constant. 
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(3) The dominant roots are a complex conjugate oair and 
the damping coefficient is to be held constant, 

(4) The dominant roots are a complex conjugate oair and 
the damped resonant or natural frequency is to be held con- 
stant. 

5.3: Case (1) : The Dominant Root is Real and is to be Held 

Constant 

When the dominant root, p^, is real and is to be held 
constant for changes in the parameters a and 3, then the 
differential of the root must be equal to zero, i.e., 



dp^ = 0. 



With 



dp^ = 0, equation 5-11 reduces to: 

, , 2 , , n 

- c-+c,p.+c~p. +...+C p. 

da _ 0 1^1 2^1 n 1 

b»+b,p.+b„p.^+. . .+b p.^ 
0 1^1 2^1 n^i 



(5-12) 



Where da/d3 is the definition of "compensating ratio". The 

constants, b. and c. are determined from the coefficients of 
3 3 

the characteristic equation. Since for this case, p^ is real, 
the compensating ratio is: 

da 



d3 



= -Constant 



(5-13) 



Example 5-1 

Figure 5-4 is a block diagram of a type one, unity feed- 
back control system. The system's characteristic equation is 
third order and contains the variable parameters a and 3. 

The characteristic equation of the system shown in figure 5-4 
is : 

s^ + (Pj^+3)s^ + (P^3+K)s + aK = 0 (5-14) 
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COMPENSATION 



PLANT 



R(s) 





(s+a) 




K 


(S+Pj^) 


s (s+8) 










Figure 5-4: Block Diagram Example 5-1 System 



Choosing the components and gain factor such that the system 
is stable let: 

Pj^ = 12; a = 3; B = 1; K = 210. 

The characteristic equation is now: 

s^ + 13s^ + 222s + 630 = 0 
which factors to : 

(s+3.318) (s+4.84+jl2.9) (s+4 . 84- j 12 . 9) (5-15) 

Without any type of compensation or self- adaption, the 
behavior of the system is dependent upon changes in the param- 
eter 8. Table 5-1 lists various values of 8 and the corres- 
ponding factors of the systems characteristic equation. For 
convenience these values are plotted in figure 5-5. 



Table 5-1: 
8 Values 
.01 
.1 
1.0 
2.0 
10.0 



System Roots With a Variable Parameter 
Factors of Characteristic Equation 
(s+3.49 ) (s+4.26+jl2.73) (s+4 . 26- jl2 . 73) 

(s+3.47 ) (s+4.31+jl2.75) (s+4 . 31- j 12 . 75) 

(s+3.318) (s+4. 84+jl2.9 ) (s+4 . 84- j 12 . 9 ) 

(s+3.1 ) (s+5.4 +jl3.0 ) (s+5.4 -jl3.0 ) 
(s+2.2 ) (s+9.9 +jl3.7 ) (s+9.9 -jl3.7 ) 
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Assuming the desired operating point for the system is 
that position specified by the value 6=1^ Equation 5-15 
gives the factors of the characteristic equation. Note that 
the dominant root of the system is = -3.318. The compensa- 
tion network has been added to provide adaption so that any 
variations in 6 can be cancelled by the adjustment of the 
parameter a insofar as the dominant root performance is con- 
cerned. 

Without specifying any change in a or 3, evaluate the 
constant term of equation 5-13, i.e., 

, . 2^ ^ n 

c^tCiPitC2Pi 

bo+bj^Pi+b2Pi^+...+b^Pi" 
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Ka 



a 



0 



bQ = K = 210; 




0 



ai = + Pi3 + K 

bi = 0; = 12; = K = 210 



b2 = 0; C 2 = 1; d2 = = 12 

as = 1 

bs = C3 = 0; ds = 1 

With the value = -3.318, equation 5-12 becomes: 
da _ 0+12(-3.318)+l(-3.318)^+0 

33 210+0 



reducing 



da _ .28.7936 
d6 210 



hence : 

da = +.13711d3 (5-18) 

From the results of equation 5-18, for any change in the 
plant's parameter 3/ in order to keep the root = -3.318 
constant, a must be changed (.137)d3. 

Situation (1) : 

Assume 3changes from a value of 1 to a new value, 3=2. 
As the change in 3 is plus one, a must be changed (.137) (1) 
to a value a' = +3.137. These new values of a and 3 change 
the characteristic equation to: 

s^ + 14s^ + 234s + 658.79 = 0 
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By a root solving method, the above equation factors to: 
(s+3.3179) (s+5. 34+jl3.04) (s+5 . 34- jl3. 04) 

The dominant root has been held constant. 

Situation (2) : 

Assume 3 changes from a value of 1 to a value 3 = 10. 

The variable parameter a must then be changed (.137) (9) = 

1.233 to a new value a' = 4.233. The new characteristic equa- 
tion becomes : 

+ 22s^ + 330s + 889.14 = 0. 

This equation factors to: 

(s+3.3174) (s+9.34+jl3.45) (s+9 . 34- jl3. 45) 

The dominant root has been held constant. 

Situation (3) : 

Assume 3 changes from a value of 1 to 3 = .9. This 10 
percent decrease forces a to change -.0137 to a new value 
a' = 2.9873. The characteristic equation becomes: 

s^ + 12. 9s^ + 220.8s + 627.123 = 0 
which factors to: 

(s+3.318) (s+4.79+jl2.88) (s+4 . 79- jl2 . 88) . 

The dominant root has been held constant. 

Example 5-2 

Given a type 'O', unity feedback system as shown in fig- 
ure 5-6. The system characteristic equation is third order 
and is given by: 

s^ + (P^+P2+3 )s^ + [P^P 2 +KK^+ (Pj^+P 2 ) 3] s + KK^a + = 0 

(5-19) 
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R(s) 





(s+a) 


, ^ 


K 




r.,(s) 


r ^ 


(s+p^) 




Ti+PjyTs+BT 

















Figure 5-6; Block Diagram Example 5-2 System. 

Let the values of the various components be assigned such that 
the system is stable. These values are: 

?! = 2; ?2 = 3; KK^ =4;a=5;6=6. 

The characteristic equation with these values substituted be- 
comes : 

s^ + lls^ + 40s + 56 = 0. 

• i 

The factors of the above equation are: 

(s+5.716) (s+2.642+jl.678) (s+2 . 642- jl . 678) . 

The results of the Case (1) compensating ratio can be 
used such that for any variation in the parameter B, one will 
be able to vary the parameter a and keep the real root , p^ = 
-5.716 constant. The procedure to be followed can be outlined 
into three steps. 

Step (1) ; Determine the value of the constant terms in the 
coefficients of the characteristic equation: 
a^ = KK^a + 

bo = = 4; Cq = P2^P2 ~ ~ ^ 



^1 " ^ 1^2 ^^1 



bi = 0; Cj^ = + ?2 = 5; d^ = Pj^P2 + = 10 



(5-20) 
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a2 = Pi + ?2 + S 






a 



3 



1 



b 



3 





(5-20) 



Step (2) : Using the values found in equation 5-20 and the 

real root = -5.716, evaluate the compensating ratio equa- 
tion 5-12. 

da _ 6+5 (-5.716)+1(-5.716)^ 

d6 Z 



da _ 
d6 “ 



9.0926 

— r~ 



or da = -2.273d6 



(5-21) 



Step (3) : Identify the change in the parameter 3 and adjust 

the value of a according to equation 5-21. The real root 
will remain constant. 

Assume 6 changes from a value of +6 to a new value of 
5.3. dB = -.7 and the change in a must be: 

da = -2.273(-.7) = 1.59 

Siabstitute the new value of a = 5 + 1.59 into the character- 
istic equation along with the new value of B. With a = 6.59 
and B = 5.3, the new characteristic equation becomes: 

s^ + 10. 3s^ + 36.5 + 58.2 = 0. 

This equation factors to: 

(s+5.69 3) (s+2.31+j2.214) (s+2 . 31- j 2 . 214 ) . 

5.4: Case (2) : The Dominant Roots are a Complex Conjugate 

Pair and are to be Held Constant 

When the dominant roots of the characteristic equation 
are a complex conjugate pair, an increment in one parameter. 
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a 



r. 



XP^ 



a 



* 



is' 



,XPi 



Figure 5-7. Change in Root Location With a 
Change in a System Parameter 



such as 3, changes the root locations as shown in figure 5-7. 

Although any new root location may occur, it is impor- 
tant to note that the two roots will always be complex con- 
jugates. The behavior of complex conjugate roots is far more 
complex than the real root case presented earlier. Instead 
of a simple plus or minus variation as in the real root loca- 
tion, an uncontrollable parameter variation in the system's 
plant will cause the complex dominant roots, p^ and to 

change both in their real and imaginary components. In this 
section, it is proven that with two variable parameters in 
a system, the dominant root locations can be returned to 
their original positions by corrective compensation in only 
very trivial cases. 

For convenience equation 5-11 is repeated below. This 
equation shows that the differential of any root, p^ , is a 
function of the change in any two variable parameters a and 

3 . 
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da + 




n 




n 



dPi 



(5-11) 



where: and Cj , j = 0,1 



n are constants determined by 



the system, 

: a and 3 are real parameters, initially fixed in 

value but subject to variations. 

If is complex and is desired to be held constant, 
then the differential of the root p^ must be set equal to 
zero , i . e . , 



As the parcimeters a and 3 are specified to be real by 
definition, then the right-hand side of equation 5-22 must 
be real to have any meaning. There must be restrictions on 
the system parameters in order to make the compensating ratio, 
equation 5-22, equal to a real constant. 

To find these restrictions, assume the numerator of the 
compensating ratio with a complex root, p^, may be combined 
into the general format: m + jn. Likewise assume that the 

denominator is expressable as: x + jy. Equation 5-21 may 

now be written: 



dp. = 0 



Equation 4-11 then reduces to: 



, c„+c, p . + . . .+c p.*^ 

da _ 0 I'^i n*^i 

d3 " n 




(5-22) 



dg _ m+jn 
d3 ~ x+jy 



(5-23) 
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Separating the above equation into its real and imaginary 
components : 



da (m+jn) (x-jv) 
38 = '(x+^y) ‘ T3RyT 



(5-24) 



or 



dg 

d8 



0 



mx+ny 

2 ^ 2 
X +y 



. nx-my 

^22 
X +y 






(5-25) 



Equation 5-25 will be real if and only if nx = my. The 



compensating ratio is : 



da 

38 



_m 

X 



_n 

'y 



(5-26) 



The numerator and the denominator of the comoensating ratio 
must be constant factors of each other. As will be shown, 
when p^ is complex the above conditions imply pole- zero can- 
cellation. This result is expressed in the following theorem. 



THEOREM 5-1 

"" th 

For the n order feedback control system shown in fig- 
ure 5-8 with variable system parameters, a and 8, there will 
exist a closed-loop root, p^ , whose location remains invar- 
ient for all increments in a or 8 such that the compensating 
ratio, da/d8, is equal to a real constant provided, if and 
only if: 

(a) the root of interest, p^ , is real or 

(b) the variable parameter a equals the variable 
parameter 8 which implies a zero-pole cancel- 
lation, a trivial case. 
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“ th 

Figure 5-8. n Order Feedback Control System 



Proof : 

For convenience, the compensating ratio for any n 
order feedback control system is repeated; 



-th 



dg 

d3 



, . . n 

Crt+c,p.+. . .+c p. 

0 I - 1 n - 1 

n 



(5-22) 



b,.+b,p.+. . ,+b p. 

0 1^1 n 1 

If the above ratio is to be equal to a real constant 'K' , 

then equation 5-22 reduces to: 

K(b-+b^p.+ +b p.^) = - (c-+c^p.+. . .+c^p.^) (5-27a) 

u 1*^1 n'^ 1 0 1 1 n 1 

equating to zero: 

(b„K+c„) + (b,K+c,)p. + + (b K+c )p.^ = 0 (5-27b) 

uu il l nn i 

A sufficient condition for part (a) of Theorem 5-1 is 
that equation 5-27 is always satisfied when the root of in- 
terest, p^, is real. This can be shown to be true by inspec- 
tion of the compensating ratio, equation 5-22. The param- 
eters a and 3 are by definition, real parameters. The con- 
stant terms, b. and c., are also defined to be real numbers. 

3 D 

If p^ is real, then equation 5-22 is a ratio of two real 
numbers, hence da/d3 is equal to a real constant. If this 
real constant is equal to 'K', then equation 5-27 must always 
be satisfied for any real root, p^. 
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A sufficient condition for part (b) of Theorem 4-1 fol- 



lows since if a = 3/ then 

|^=K=1 (5-28) 

The closed loop system in figure 5-8 is independent of a and 

3. 

The necessary condition for the proof of Theorem 5-1 can 
be shown with the aid of figure 5-9 where the real factor, k, 
has been added for generality. 



Figure 5-9. Modified Closed Loop Feedback Control System 

The closed-loop control system of figure 5-8 has been 
transformed into the modified control system of figure 5-9 
by the following modifications: 

Gj^(s) = G(s)/H(s), the fixed portion of the closed 
loop system. 

G 2 (s) = k (s+a) / (s+3) , the adjustable portion of the 
closed loop system. 

The characteristic equation is; 





> ..ciaj 




(5-29a) 



or 



H(s) (s+3) + kG(s) (s+a) = 0 



(5-29b) 



or 



s [kG (s) +H (s) ] + akG(s) + 3H(s) = 0 



(5-29c) 
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(5-30a) 



Assume that G(s) and H(s) may be expressed as: 

G(s) = go + giS + ... + gn-i®””^ 

and 

H(s) = hrt + h^s + ... + h ,s^“^ (5-30b) 

U i n-i 

Substituting equations 5-30a and 5-30b into equation 5-29c, 
the characteristic equation becomes: 

(kg^a+h^B) + (kg^a+h^^B+kg^+hQ) s + ... 

+ (kg T ot+h ^ 3+kg ^+h s 
^n-1 n-1 ^n-2 n-2 

+ (kg ,+h Js" (5-31) 

^n-1 n-1 

““ th 

Equation 5-31 is the n order characteristic equation of 
the feedback control system of figure 5-8, which is written 
as : 

a,. + a,s + ... + a , s^ ^ + a s^ = 0 (5-32) 

0 1 n-1 n 

From the introduction, the coefficients of the characteristic 

equation are expressable as: 

a. = b.a + c.B + d. (5-4) 

3 D D 3 

j = 0,1,. . . ,n 

Equating the coefficients of equations 5-31 and 5-32 yields: 
a^ = kg^a + h^B = b^a + c^B + d^ (5-32a) 

^1 ~ kg^a + h^B + (kg^+h^) = b^a + c^^B + d^^ (5-32b) 

• • • 

a ,=kg ,a + h i6 + (kg _)=b ^a + c ^B + d , 

n-1 ^n-1 n-1 ^n-2 n-2 n-1 n-1 n-1 

(5-32c) 

a„ = (kg ,+h ^) = b a + c B + d (5-32d) 

n n-i n-1 n n n 
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Equating the constant terms: 



= kg. 



= h. 



= 0 (5-34a) 



= kg^ 



Cl = hi 



di = kq^ + hg (5-34b) 



b 



n=l ^^n-] 



n 



= 0 



1 — 1 

1 

c 

II 

1 — 1 
c 


d T 

n-1 


II 

1 

to 


h . 
n-2 


n = ° 


d 

n 


= kg . + 
^n-1 


h T 

n-1 



With the above values of b^ and c^ , equation 5-27b becomes: 
(kg^K+hJ + (kgiK+hi)pi + ... + (kgj^_iK+h^_i) = 0 (5-35) 





G(s) 

HliT 




K 



















+• Vi 

Figure' 5-10. (n-1)^ Order Feedback Control System 



Figure 5-10 shows an (n-1)^ order closed loop feedback 
control system. The characteristic equation of this system 
is : 

Irlf 1 

or 

KG(s) + H(s) = 0 (5-36b) 

If G(s) and H(s) are represented by equations 5-30, then 
equation 5-36b can be written as: 

(g^K+hQ) + (giK+hi)s + ... + ^ ° (5-37) 
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Comparing equation 5-35 with equation 5-37, it follows that 
the root, p^, must also be a root of the system shown in fig- 
ure 5-10, that is, 

kK ||||- +1 = 0 (5-38) 



Since p^ is a root of the original system shown in figure 
5-8, then it also follows that; 



(p.+a) V 

' G (s) 

TpT+FT hTsT 



+ 1 



0 



Equating equations 5-38 and 5-39: 



(p^+a) 
(Pi+3) " 



Solving for the value of the root, p^^ : 

(p^+a) = kK(p^+3) 



(5-39) 



(5-40) 



(5-41) 



or 



6-a 

1-kK 



(5-42) 



From equation 5-42, the root of interest, , is always 
real if kK is not equal to one. If kK is equal to one, then 
p^ is indeterminate, however from equation 5-28, kK = 1 im- 
plies a = B or zero-pole cancellation, a trivial case of 
compensation . 

5.5: Case (3) : The Dominant Roots are Complex Conjugates and 

the Damping Coefficient is to be Held Constant 
Since the damping coefficient, is to be held con- 

stant for variations in the plant's parameter 3, let the real 
part of either dominant complex conjugate root be held con- 
stant , i . e. , 




CONSTANT 



(5-43) 
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With equal to a constant value the real part of 

the differential of must equal zero, or 

(5^3 [dp.] = 0 (5-44) 

With the above condition, the real part of equation 5-11 is 
written : 



■ • -Ig.[ 



da 



d3 



(5-45) 



hence : 






bo+biPi+.--+Vi 

-F' ip^l 



n 



da 



h-e 



CQ+CiPi+...+C^Pi 

_ -F' (P^) 



de 



(5-46) 

Evaluating the terms in the parenthesis of equation 5-46, let 



m + jn = 



and 



u + ]V = 



+ * • *+b P- 
-F' (p^) 



Co+CiPi+. . *+Cj^p^ 

-F' (p.) 



n 



n 



(5-47) 



(5-48) 



Using the above notation, equation 5-16 is rewritten as: 

(J^^[mdc.+jndol = [ud3+jvd6] (5-49a) 



or 



mda = ud3 (5-49b) 
The compensating ratio for a constant damping coefficient is 
therefore equal to: 



da _ u 
d3 m 



(5-50) 
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Example 5-3 



Using the system of example 5-2, the characteristic egua- 
tion is: 

S^ + (Pj^+P2 + B)s^ + [P^P2+KK^+ (P^+P2) 3] S + KKj^a + ° 

with the parameter values: 

Pi = 2; P 2 = 3; KK^ =4;a=5;3=6. 

The numerical characteristic equation is 

s^ + lls^ + 40s + 56 = 0 
which factors to; 

(s+5.716) (s+2.642 + jl.678) (s+2 . 642- jl. 678) . 

The parameter 3 will vary and again must be compensated 
for by a change in the parameter a, such that the damping 
coefficient, is held constant. The damping coefficient 

will be held constant if equation 5-50 is satisfied, there- 
fore the following procedure is recommended. 

Step (1) ; Determine the value of the constant terms in the 
coefficients of the characteristic equation. 

This step has been performed in example 5-2 and equa- 
tion 5-20 is repeated for convenience. 

bo = 4 Co - 6 <So = ° 

bj^ = 0 Cj^ = 5 dj^ = 10 

b2 = 0 C 2 = 1 ^2 = ^ 

bs = 0 C3 = 0 ^3 = 1 

Step (2) : Form the real part of equations 5-47 and 5-48. 
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Denominator : 



3 2 

F(s) = s + 11s + 40s + 56 

F' (s) = 3s^ + 22s + 40 

= -2.64 - jlO.35 

F' (p^) = 3(-2.64-jl.68)^ + 22 (-2 . 64-jl . 68) + 40 
F' (p^) = -5.64 - jlO.35 



Equation 5-47; 



m + 



in = 



-F' (Pi> 



m = 



5.64+]10.35 

4(5.64) 

(5.64)^+(10.35)^ 



Equation 5-48: 



m = .163 



C +C p.+C^p. 

“ ^ (p. ) 

_ -3.1+j.43 _ 3.13 /7 . 9 

- 5.64+jl0.35 - 11.8 /61.T 

u = -.094 



The compensating ratio is the negative ratio of u over m or 

da = +.573d8 

Step (3) : Identify the change in the parameter 8 and adjust. 

Assume 6 changes from a value of 6 to a value of 5.3. 

The change in 6 is; 

d3 = -.7 

da = (-.7) (.573) 

da = -.4011 
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The parameter a must change to 5 - .4011 = 4.598^, The new 
characteristic equation will be: 

+ 10. 3s^ + 36.5s + 50.2 = 0. 

This equation factors to: 

(s+5.01) (s+2.64+jl. 74) (s+2 . 64- j 1 . 74) . 

By inspection, the real part of the complex conjugate pair 
has been kept constant. 

5.6: Case (4): The Dominant Roots are a Complex Conjugate 

Pair and the Damped Resonant Frequency is 
to be Held Constant 

In a similar manner as the damping coefficient, was 

held constant in case (3), the damped resonant or natural 
frequence, will be held constant by keeping the imaginary 

part of either dominant complex root constant, i.e., 

^^mtPiJ = CONSTANT 



With the Imaginary part of p^^ held constant, the imaginary 
part of the differential of p^ must equal zero, or 



^„(dPi] = 0 (5-51) 

Setting the imaginary part of the differential of p^ equal 
to zero, equation 5-11 becomes: 






m 



-F' (p^) 



n 



da 



■ c^m 



_ -F' (p^) 



n 



• d6 
(5-52) 



Inspecting the terms in the parenthesis as done in case (3), 
equations 5-47 and 5-48 are evaluated. The comoensatina 
ratio for a constant damped resonant frequency is found to be 

(5-53) 



da 

3F 



V 

n 
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where 



V 






Cr, + C,p. + . , .+C p. 
0 I'^i n^i 

-F' (p^) 



n 



(5-54) 



and 



n = 



b^+b-p,+. . ,+b p. 
0 I'^i n^i 



n 



-F’ (p.) 



(5-55) 



Example 5-4 

Using the system of example 5-2 and 5-3, the characteris- 
tic equation is : 

s^ + (Pj^+P2+6)s^ + [Pj^P2+KKj^+(P^+P 2) els + KK^^a + Pj_P2^ " ° 

Substituting in the values of the parameters as before, the 
characteristic equation becomes: 

s^ + lls^ + 40s + 56 = 0 
which factors to: 

(s+5.716) (s+2.64+jl.68) (s+2 . 64- jl . 68) . 

In this example the damped resonant frequency will be 
held constant by the following procedure as indicated in the 
above case. 

Step (1) : Determine the value of the constant terms in the 

coefficients of the characteristic equation. 

Repeating equation 5-20 for convenience: 



bo = 4 



Co = 6 



d,, = 0 

0 



bj^ = 0 



Cj^ = 5 



= 10 



b2 = 0 



^2 ~ ^ 



d2 = 5 



b3 = 0 



C3 = 0 



d3= 1 



Step (2) : Fom the imaginary parts of equation 5-47. 
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Denominator : 



F(s) = + lls^ + 40s + 56 

F' (s) = 3s^ + 22s + 40 

= -2 . 64- jlO . 35 : the dominant root of 
interest. As before in example 5-2: 

F' (p^) = -5.64-jl0.35. 

Equation 5-47; 

4 

“ 5. 64+jlO. 35 
4 

“ 11. 8 /61. 5~ 
n = - . 298 

Equation 5-48; 

_ -3.1+j.43 
5.^4+jl0.35 

= .266 /110°46' 

V = .246 

Step (3) ; Determine the compensating ratio according to 
equation 5-53. 

da _ .246 

d6 " ”-.298 

da = +.825dB 

Step (4) : Identify the change in the parameter 3 and adjust. 

If 3 changes from a value of 6 to a new value 3 = 5.3, 
the change in 3 is: 

d3 = -.7 

da = (-.7) (.825) 

da = -.577 
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The parameter a must change to 

5 - .577 = 4.423 . 

The new characteristic equation will be: 

s^ + 10. 3s^ + 36.5s + 49.4 = 0. 

This equation factors to: 

(s+4.94) (s+2.68+jl.683) (s+2 . 68- jl. 683) . 

By inspection the imaginary part of the dominant root 
has been held constant by proper use of the compensating ra 
tio, equation 5-53. 
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CHAPTER VI 



Conclusions 

In this thesis, a purely mathematical property of all 
polynomials of degree n, (called the associated oolynomials 
of the defining polynomial) has been applied to the fields 
of circuits and systems. The properties of associated polyno- 
mials have been investigated and examples given that demon- 
strate that these polynomials provide useful results. It 
has been shown that associated polynomials may be used in a 
new method for accomplishing a partial fraction expansion in- 
cluding the troublesome case of repeated roots of high order 
of multiplicity. Associated polynomial theory has been used 
to study the relationship between variations of polynomial 
coefficients and polynomial roots. The theorems developed 
in Chapter Three that show that the variation of the coeffi- 
cients of a polynomial with respect to the variation of its 
roots or vice-versa is determined by associated polynomials 
have important engineering applications. They are immedi- 
ately used in a root solving procedure and to express sensi- 
tivity coefficients in a new analytical form. 

The most significant result of this thesis is the 
investigation of the theory of compensating parameters as 
applied to self-adaptive control systems. A mathematical 
relationship between two variable parameters^ such -diat the 
system's dominant root performance, has been derived. In 
addition, a theorem of fundamental importance in adaptive 



91 



control theory has been presented that Effects the basic 
theory of complex- conjugate dominant root adaptive compensa- 
tion. 

As a suggestion for further research, it is recommended 
that an actual computer program be written that performs the 
partial fraction expansion of any ratio of polynomials. In 
addition, it is believed that the results of Chapter Three 
have logical extensions to a derivation of an analytical ex- 
pression for a root-locus plot. This analytical expression 
could then be programmed onto a digital computer. Lastly, 
the results of Chapter Five can be extended to variable 
parameters that are not linearly related to the coefficients 
of the characteristic equation, i.e., the variable parameters 
would appear in the form: 

a. = b.a + C.3 + h.a6 + d. (6-1) 

3 3 3 3 3 

or : 

2 2 
a. = b.„a + b.,a + h.a3 + c.,B + c.„B + d. (6-2) 

T t2 tI t t1 t2 t 
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APPENDIX A 



Horner's Method 

A.l: Introduction 

A complete historical survey of associated polynomials 

is given in reference 1. In this reference it is noted that 

what has been labeled associated polynomials are sometimes 

referred to as Horner polynomials. This appellation comes 

about by three processes involving associated polynomials 

(4) 

that are attributed to Horner . These processes are: 

(1) : The nested evaluation of polynomials; 

(2) : Synthetic division; 

(3) : A method of approximating the roots of polyno- 

mial equations. 

Each of the three processes above are commonly called Horner's 
method or rule^^^^. As is shown, the method of approximating 
the roots of a polynomial directly involves synthetic divi- 
sion. 

A. 2: The Nested Evaluation of Polynomials 

Associated polynomials may be used in the evaluation of 
polynomials by a process called "nesting". In order to eval- 
uate the polynomial, 

F(s) = a s^ + a - s^ ^ + ... + a,.: a =1 (1-1) 

n n-1 0 n 

for any value of s , there is required 2n + 1 multiplications 
and n additions. Using associated polynomials there is re- 
quired only n multiplications and n additions. This process 
of evaluation is performed by the following method: 

F(s=u) = a u^ + a iU^"^ + . . . + a^ : a =1 (A-1) 

n n-1 u n 
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but 



F(u) = Aq(u) 

where ; 

Aq(u) = uA^(u) + 

Aj^(u) = uA2(u) + 



( 1 - 12 ) 

(l-9n) 

(l-9n-l) 



and since: 



A , (u) = uA^ (u) 
n-1 n 



+ a, 



n-1 



(l-9b) 



A„^^(u) = 0 (1-7) 

A^(u) = nA^+i<n) 

Hence F(u) is evaluated using only n multiplications and n 
additions . 

A. 3: Synthetic Division 
The equation: 



F(s) 

s-u 



A (u) 



I A. (u)s^ ^ 

j = l ^ 



(1-15) 



states that the quotient and remainder after division by a 
polynomial of the form s-u may be calculated recursively. 

The nature of the recursion is such that one may work with 
detached coefficients. By detached coefficients is meant an 
abbreviation of the ordinary multiplication and division 
processes used in ordinary algebra. The coefficients alone 
with their signs are used in the process. The powers of the 
variable occurring in the various terms are understood only 
from the order in which the coefficients are written. Missing 
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powers are assumed to be present with zero coefficients. The 
process of synthetic division is best understood by the follow- 
ing example and summary. 

Example A-1 

Given the polynomial; 

F(s) = As^ + BS^ + Cs + D (A-2) 

And the divisor polynomial: 

s - a (A-3) 

The usual process of division is as follows: 

As^ + Bs^ + Cs + D / s-a 

As^ - Aas^ As^+Es+F 

2 

Es + Cs 
2 

Es - Eas 

Fs + D 
Fs - Fa 

D + Fa = R 



Or ; 



3 2 

As +Bs +Cs+D 2.„ . 

= (As +Es+F) + 

S "* 3, 



S-a 



(A-4) 



Comparing with equation 1-15; 

AQ(a) = R (A-5) 

and 

Aj(a)s^ ^ = (As^+Es+F) : j = l,2,...,n (A-6) 

Using synthetic division, the process of division is reduced 
by only using the coefficients according to the following 
scheme . 



A B C D /+a 

Aa Ea Fa 

A E F R 
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Notice that the first three mimbers of the last line, 

2 

A, E, and F, are the coefficients of s , s, and 1, in the auo- 

tient, and R is the remainder. The process may be exoanded 
th 

to any n order polynomial. 

As a summary, synthetic division may be performed accord- 
ing to the following rule. 

(1) : To divide F(s) by a polynomial of the form, s - a, 

arrange F(s) into descending powers of s. 

(2) : Write the coefficients of F(s) on a horizontal 

line, in the order which corresponds to the arrangement spec- 
ified in (1). If any power of s is missing in F(s) , supply 
that power with a zero coefficient. 

(3) : Multiply the first coefficient A by a, write the 

product below the second coefficient B and add. Multiply 
this sum, E, by a, write the product below the third coeffi- 
cient, D, and add. 

(4) ; Proceed in this way until all of the places in the 
third row except the first are filled up, and write down the 
first coefficient. A, of F(s) in the first place of the third 
row. The last number of the third row will be the remainder, 
and the other numbers in the third row will be the coeffi- 
cients of the quotient obtained when F(s) is divided by s - a. 
A. 4: Approximation of the Roots of a Polynomial Equation 

A method for approximating the real roots of an algebraic 

equation first appeared in a paper read by W. G. Horner be- 

( 4 ) 

fore the Royal Society in 1819 . This process involves 

associated polynomials in that repeated use of synthetic 
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division is used to expedite the work of substituting values 

( 21 ) 

for the variable. The essential steps of Horner's method 

of approximation are as follows: 

(1) : Using curve plotting, the theorem of continuous 

fractions*/ etc., isolate a positive root between two succes- 
sive integers. If the equation has only negative roots, 
transform it to one whose roots are the negatives of those 

of the given equation. 

(2) : Transform the equation into an equation whose roots 

are decreased by the lesser of the integers between which the 
root lies, by the substitution x' = x-a. The root of the new 
equation will lie between zero and unity. 

(3) : Isolate the root of the new equation between 

successive tenths. 

(4) : Transform the last equation into an equation whose 

roots are decreased by the smaller of these tenths and iso- 
late the root of this equation between hundredths. Continue 
this process to one decimal place more than the place to 
which the answer is to be correct. 

The root sought is then the total amount by which the 
roots of the original equation were reduced, namely, the 
lesser integer, plus the lesser tenth, olus the lesser hun- 
dredth, etc., the last decimal being rounded off to make the 
result accurate to the desired decimal place. 

*If P(x) is a polynomial with real coefficients and of 
degree n>l, and if a and b are real numbers such that P(a) 
and P(b) have opposite signs, then the polynomial P(x) has 
at least one root between a and b. (22) 
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Example A-2 



Find to two decimal places that root of the equation: 



s^ + - 12s + 11 = 0 (A-7) 

which lies between 2 and 3. 

The root in question will be of the form 2.abc..., 
where a, b, c, ... stand for its decimal figures. 

(1) : Using synthetic division, first transform equation A-7 

into an equation whose roots are those of equation A-7 dimin- 
ished by 2. 



1 1 -12 11 /2 

6 -12 
13 - 6-1 
2 10 

I 5 T 

2 

1 7 



The transformed equation is: 

s^ + 7s^ + 4s - 1 = 0 (A-8) 

Equation A-8 will have the root .abc... lying between 0 and 

1 . 

(2): Sxibstituting the values s= .1, .2, .3, ... into equa- 

tion A-8, one finds that: 

f(.l) < 0 

and 

f(.2) > 0 

hence equation A-8 has a root between .1 and .2, so that: 

a = 1 (A-9) 
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(3): Diminishing the roots of equation A-8 by .1 using syn- 

thetic division: 



17 4 -1 /.I 

.1 .71 .471 

I TTI 4T71 - .529 

.1 .72 

I 772 5TT3 

.1 

1 7.3 



The transformed equation is: 

s^ + 7.3s^ + 5.43s - .529 = 0 (A-10) 

which has a root of the form . Obc . . . 

(4): Substituting the values s = .01, .02, ... into equa- 

tion A-10, one finds that: 

g ( . 08) < 0 

and 



g(.09) > 0 . 

Hence equation A-10 has a root between .08 and .09 so that: 

b = 8 . 



From equation A-9 and A-11, the root of equation A-7 
which lies between 2 and 3 to two decimal places is : 

s = 2.18 . 

The process could be continued indefinitely, thus per- 
mitting one to compute the root of equation A-7 between 2 and 
3 to any specified number of decimal places. 
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